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'This  report  is  a revised  and  extended  version  of  a previous  study  entitled  "On 
the  Theory  of  Steady  Motion  of  Low-Froude-Number  Displacement  Ships"  that  was 
submitted  and  accepted  for  publication  in  the  Journal  of  Ship  Research  in 
October  1977  and  March  1978,  respectively.  Publication  of  this  study  however 
was  interrupted  due  to  its  reliance  upon  a preliminary  study  entitled  "Potential 
Theory  of  Steady  Motion  of  Ships"  that  is  still  in  the  process  of  being  reviewed 
for  publication.  An  updated  version  of  the  latter  study  has  been  presented  in 
two  previous  reports  [1,2],  which  — together  with  the  present  report  — provide 
a coherent  and  detailed  exposition  of  the  "generalized  Neumann-Kelvin  theory" 
of  ship  wave  resistance.  \ , jit 


ABSTRACT 


The  main  object  of  the  present  Part  4 of  the  "generalized  Neumann-Kelvin 
theory"  of  ship  wave  resistance  is  to  present  "low-Froude-number  wave-resistance 
approximations"  as  well  as  "low-Froude-number  slender-ship  approximations",  which 
supplement  the  basic  "slender-ship  approximations"  obtained  in  Parts  2 and  3. 
First-  and  second-order  approximations  are  given  explicitly.  In  particular, 
the  low-Froude-number  wave-resistance  approximations  of  Guevel  et.  al. , Kayo, 
Baba,  and  Maruo  are  obtained  in  this  study  as  first-order  approximations  in  a 
sequence  of  low-Froude-number  iterative  approximations  for  solving  the  basic 
"generalized  Neumann-Kelvin  integral  equation"  established  in  Part  2 of  the 
theory.  New  "first-  and  second-order  low-Froude-number  slender-ship  wave- 
resistance  approximations"  are  derived.  In  addition,  a new  integral  equation 
for  evaluating  potential  flow  about  a ship  hull  in  the  zero-Froude-number  limit" 
is  obtained.  A method  of  solution  of  this  integral  equation,  based  on  an 
iterative  procedure  leading  to  a sequence  of  "slender-ship  approximations",  is 
discussed  in  some  detail. 
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INTRODUCTION 

In  this  study,  the  exposition  of  the  "generalized  Neumann-Kelvin  theory"  of 
ship  wave  resistance  developed  in  [1]  and  [2]  is  pursued.  Specifically,  the  object 
of  the  present  Part  4 of  the  theory  is  to  derive  "low-Froude-number  approximations", 

• as  well  as  "slender-ship  low-Froude-number  approximations",  which  supplement  (and  in 

a sense  generalize)  the  basic  "slender-ship  iterative  approximations  given  in  [1]  and 
[2] . The  plan  and  main  results  of  the  study  are  briefly  described  below. 

In  the  first  section,  the  basic  differential  equations  of  the  problem  are  listed, 
and  the  fundamental  integral  equation  underlying  the  "generalized  Neumann-Kelvin  theory" 
developed  in  [1]  and  [2]  is  recalled.  Section  two  is  concerned  with  the  analysis, 
and  more  particularly  with  the  numerical  evaluation,  of  the  velocity  potential,  4>  » 
of  the  disturbance  flow  caused  by  the  ship  in  the  "zero-Froude-number  limit".  The 
well-known  result  that  the  "zero-Froude-number  limit",  0 , of  the  disturbance  potential 
4>  is  identical  to  the  "rigid-wall  potential"  — in  which  the  free  surface  is  replaced 
by  a horizontal  rigid  wall  coinciding  with  the  undisturbed  free  surface  — is  first 
established  (in  two  different  ways,  by  starting  from  the  differential  and  integral 
formulation  of  the  problem).  The  main  object  of  section  two,  however,  is  to  present 
' a new  integral  equation,  namely  equation  (17) , for  the  numerical  evaluation  of  the 

"zero-Froude-number  potential"  Two  alternative  equivalent  forms,  given  by  equa- 

tions (19)  and  (25) , of  a recurrence  relation  for  solving  the  "zero-Froude-number 
integral  equation"  (17)  iteratively  are  given.  These  equivalent  recurrence 
formulas  generate  a sequence  of  "slender-ship  zero-Froude-number  approximations" 
tf’k*  k ^ 0,  where  the  "zeroth  approximation"  is  taken  as  zero,  i.e.  4>q  =0,  which 
becomes  exact  in  the  limit  of  an  extremely  slender  ship-hull  form,  and  the  "first 
approximation"  <p®  is  the  potential  defined  by  formula  (18)  and  referred  to  as  the 
"initial  approximation".  This  "initial  approximation"  is  examined  in  some  detail. 

In  particular,  it  is  shown  that  the  approximation  4)^  can  be  regarded  as  a generalization 
of  the  classical  "thin-ship  approximation",  to  which  it  reduces  in  the  limit  of  a thin 
ship.  It  is  also  shown  that  in  the  case  of  a ship  hull  in  the  form  of  a half  ellipsoid 
(say  with  beam/length  and  draft/length  ratios  equal  to  6 and  6,  respectively)  the 
. Initial  potential  <J>^  actually  is  proportional  to  the  exact  zero-Froude-number  potential 

<p0,  that  is  we  have  <p^  = k(B,<5)<j>!j!  where  the  multiplicative  correction  function  k(8,<$) 

» is  defined  by  formula  (28a)  and  depicted  in  figure  2.  For  8 - .2  and  6 * .1,  which 

can  be  regarded  as  typical  values  of  the  beam/length  and  draft/length  ratios,  we  have 
k(8,6)  1.06,  so  that  the  initial  potential  <J>^  given  by  formula  (18)  in  fact  provides 

a fairly  accurate  approximation  to  the  zero-Froude-number  potential  of  the  disturbance 
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flow  due  to  an  ellipsoidal  hull  form  of  typical  ship-like  dimensions.  To  be  sure, 

real  ship  hulls  are  not  ellipsoidal  in  shape.  Nonetheless,  the  initial  potential 

[the  modified  potential  k(3,6)^  may  also  be  used]  seems  likely  to  provide  a realistic 

approximation  to  the  zero-Froude-number  disturbance  potential  <$P . 

In  section  three,  "first-order  low-Froude-number  wave-resistance  formulas"  are 

obtained  by  approximating  the  near-field  disturbance  potential  <j>  in  expression  (31) 

for  the  Kochin  free-wave  spectrum  function  ft(t)  by  the  "zero-Froude-number  potential" 

$0.  The  basic  expression  for  the  "first-order  low-Froude-number  approximation", 

2F 

ft.  (t),  to  the  Kochin  spectrum  function  ft(t)  is  given  by  formula  (37).  An  alternative 

ZF 

expression  for  the  first  approximation  ft^  is  given  by  formula  (40) , which  has  been  ob- 
tained from  formula  (37)  by  making  use  of  a classical  Green  identity.  This  alternative 

2.F 

expression  for  the  first  approximation  then  readily  yields  the  low-Froude-number 
approximations  obtained  previously  by  Guevel,  Vaussy,  and  Kobus  [3],  Kayo  [4],  Baba  [5], 
and  Maruo  [6],  which  thus  are  obtained  in  this  study  as  first-order  approximations  in 
a sequence  of  low-Froude-number  iterative  approximations  (the  second-order  approxi- 
mation in  this  sequence  is  presented  in  section  five)  for  solving  the  basic  "gener- 
alized Neumann -Kelvin  integral  equation"  (6).  If  the  slender-ship  approximation 
k(B,6)<j)j  to  the  zero-Froude-number  potential  $ is  used,  instead  of  $ , as  an  approxi- 
mation to  the  near-field  disturbance  potential  <J>  in  expression  (31)  for  the  Kochin 

spectrum  function  ft(t) , the  "first-order  low-Froude-number  slender-ship  approximation" 
JtFsi 

ft^  may  be  obtained.  This  approximation  is  expressed  by  formula  (45)  or  equivalently 
by  formulas  (46)  and  (47).  The  simplest  (and  thus  perhaps  the  most  attractive  for 
practical  purposes)  of  the  "first-order  low-Froude-number  approximations"  obtained  in 

section  three  is  the  "linearized  (that  is,  the  nonlinear  free-surface  correction  terms 

ilFsJl 

are  neglected)  first-order  low-Froude-number  slender-ship  approximation",  ft^  , 

given  by  formula  (48).  In  fact,  comparison  of  this  formula  with  formula  (32)  for  the 

basic  "zeroth-order  slender-ship  approximation"  ft.  shows  that  formula  (48)  for  the 
iFsi  U 

approximation  ft^  is  not  significantly  more  complex  than  formula  (32)  for  ftg  from 
the  computational  point  of  view. 

In  section  four,  various  (closely-related)  "first-order  low-Froude-number  approxi- 

mations"  <j>1  and  corresponding  "first  low-Froude-number  slender-ship  approximations" 

<t>i  to  the  disturbance  velocity  potential  4>  are  obtained  by  approximating  the  unknown 

potential  4>  in  the  "Neumann-Kelvin  integral  equation"  (6)  by  the  "zero-Froude-number 

potential"  iP  and  the  corresponding  slender-ship  approximation  k(B,<$)4>j»  respectively. 

This  section  thus  parallels  and  supplements  section  three;  specifically,  the  "first- 

order  low-Froude-number  approximations"  <{),  and  4>.  to  the  disturbance  velocity 

11  JlF 

potential  $ correspond  to  the  "first-order  low-Froude-number  approximations"  ft^  and 
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£2..  to  the  Kochtn  free-wave  spectrum  function  £2  obtained  in  section  three,  and 
^ J,p  jj,FgJ2, 

indeed  the  approximations  £2  and  £2.  could  alternatively  have  been  derived  from 

^ jlFs£ 

the  analysis  of  the  "far-field  limit"  of  the  approximations  <j>^  and  to  the 

disturbance  potential. 

£p  j^p 

In  section  five,  "second-order  low-Froude-number  approximations"  and  £2„ 

JlFsJC  JlFsil 

and  corresponding  "low-Froude-number  slender-ship  approximations"  ^ and  £22 

to  the  disturbance  velocity  potential  <}>  and  the  Kochin  free-wave  spectrum  function 

£2  are  obtained  by  using  the  above-defined  "first-order  low-Froude-number  approxi- 
£p  J»Fsil 

mations"  and  (J)^  as  approximations  to  the  potential  <}>  in  the  "Neumann-Kelvin 

integral  equation"  (6) . Of  greatest  interest  for  practical  purposes  may  be  the 

£FsJl 

"second-order  low-Froude-number  slender-ship  wave  spectrum  approximation"  £22 
defined  by  formula  (64) , which  is  based  upon  the  "linearized  first-order  low-Froude- 


number  slender-ship  potential"  4>i 


given  by  formula  (58) . While  the  second  approxi- 


mation £2„ 


evidently  is  considerably  more  complex  than  the  first  approximation 


£ps£ 

£2^  given  by  formula  (48) , it  nonetheless  appears  to  be  well  within  present-day 
calculation  capabilities.  Indeed,  the  major  computational  task  is  that  associated 

£Fg  £ 

with  the  evaluation  of  the  potential  defined  by  formula  (58) , which  is  not 

significantly  more  complex  than  the  "initial  potential"  4)^  given  by  formula  (7)  or 

than  the  classical  Michell  potential.  Naturally,  a significant  simplification  of 

£Fs  a 

formula  (64)  for  the  approximation  £2»  can  be  achieved  if  the  hull-flux  correction 
£.Fs£.  ^ 

term  q,  , 0 , associated  with  sinkage  and  trim  effects,  and  the  nonlinear  free-surface 

tllJt  i-pgSL 

correction  term  are  neglected. 

In  section  six,  it  is  briefly  indicated  how  the  low-Froude-number  approximations 
4>^F,  and  ^ — which  have  been  obtained  in  the  previous  sections  as 

first  and  second  approximations  in  a sequence  of  iterative  approximations  for  solving 
the  basic  "generalized  Neumann-Kelvin  integral  equation"  (6)  — can  also  be  obtained 
by  using  two  alternative  (but  equivalent  and  closely-related)  approaches. 
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1.  The  generalized  Neumann-Kelvin  theory 

The  problem  examined  in  this  study  is  that  of  flow  caused  by  a displacement  ship 
in  steady,  rectilinear  motion  at  the  free  surface  of  an  otherwise  calm  sea,  which  is 
assumed  to  be  of  infinite  depth  and  lateral  extent.  Water  is  supposed  to  be  homogeneous 
and  incompressible.  Surface  tension  is  neglected.  Viscosity  effects  are  ignored  and 
irrotational  flow  is  assumed.  Effects  of  wavebreaking  and  spray  formation  at  the  ship 
bow  are  also  neglected.  (The  ad-hoc  corrections  for  effects  of  viscosity,  spray,  and 
wavebreaking  which  were  included  in  the  potential-flow  theory  developed  in  [1,2]  will 
not  be  incorporated  in  this  study  for  shortness.) 

The  flow  is  observed  from  a moving  system  of  coordinates  attached  to  the  ship,  so 
that  the  flow  becomes  independent  of  time.  The  z axis  is  chosen  vertical,  positive  up- 
wards, with  the  undisturbed  free  surface  taken  as  the  plane  z = 0.  The  x axis  is  par- 
allel to  the  direction  of  motion  of  the  ship  and  positive  towards  the  ship  stem.  Flow 
variables  are  rendered  dimensionless  with  respect  to  some  reference  length  L,  the  speed 

of  the  ship  U,  and  the  density  of  water  p.  The  reference  length  L may  be  chosen  as 

2 

L = U /g,  where  g is  the  acceleration  of  gravity,  as  in  [1,2],  although  it  may  be  more 
convenient  to  take  L as  the  length  of  the  ship  in  the  present  study,  which  is  specif- 
ically concerned  with  the  steady  motion  of  a ship  at  low  Froude  number.  The  parameter 


U/(gL) 


1/2 


is  denoted  by  the  symbol  F,  which  represents  the  usual  Froude  numbar  if  L is 


the  length  of  the  ship,  while  we  simply  have  F = 1 if  the  reference  length  L is  taken 

2 -V 

as  U /g  as  in  [1,2].  We  thus  define  the  dimensionless  coordinates  x = X/L  and  velocity 
potential  4>  = 4>/UL,  where  X and  $ are  dimensional,  while  the  dimensionless  velocity  is 
given  by  <4>x,4>y ,4>?)  = (*X,*Y,*Z> /U. 

The  hydrodynamic al  problem  amounts  to  determining  the  (dimensionless)  velocity 
potential  <b  of  the  disturbance  flow  caused  by  the  ship.  As  it  is  explained  in  [1], 
this  problem  may  be  formulated  in  a "solution  domain",  (d)  say,  bounded  by  some  arbi- 
trary "fictitious  hull"  surface,  (h)  say,  which  may  for  instance  (but  need  not)  be 
taken  as  the  wetted  hull  of  the  ship  in  position  of  rest,  and  the  "undisturbed  free 
surface",  (f)  say,  defined  as  the  portion  of  the  plane  z = 0 located  outside  the  inter- 
section curve  (c)  of  (h)  with  the  plane  z * 0.  The  problem  then  consists  in  solving 
the  Laplace  equation 


7 <t>  - 0 in  (d) 


(1) 


subject  to  the  "Kelvin  boundary  condition" 


IB 


5 


<f>2  + F ^ - - qf  on  (f)  , (2) 

the  "radiation  condition"  of  no  waves  upstream  from  the  ship,  and  the  "Neumann 
boundary  condition" 

<j>n  * - V - qh  on  (h)  , (3) 

where  <f>n  = V(f>*n  and  V = n • i,  with  n representing  the  unit  normal  vector  to  (h)  di- 
rected  towards  the  interior  of  the  ship,  and  i defined  as  the  unit  positive  vector 
along  the  x axis. 

The  "f ree-surface  flux"  in  the  free-surface  boundary  condition  (2)  accounts 

for  both  the  nonlinear  terms  in  the  free-surface  condition  and  the  fact  that  the 

2 2 

position  of  the  free  surface,  defined  by  z = -F  (4>x  + |V4>|  /2)  , generally  differs  from 
the  plane  z = 0 (the  undisturbed  free  surface)  where  the  free-surface  condition  is 
enforced  for  mathematical  simplicity.  The  nonlinear  free-surface  correction  flux 
q^  is  given  by 


qf  - <PZ  + F2^  + ( |V4>|  2)x  + |V4>*V|  V<b 1 2>J 


z = -Y  (<J>x  + y!V<t>|  ) 


1 9 “ <t> 

+ 4lvd)|2)  L 


4)  + F <t> . 


The  "hull  flux"  q^  in  the  hull  boundary  condition  (3)  is  associated  with  the 
fact  that  the  fictitious  hull  surface  (h)  where  the  hull  condition  is  enforced  may 
differ  from  the  actual  ship  hull  surface,  (H)  say.  The  hull  flux  q^  is  given  by 

= (?  + V<J>)H  • N - (i  + 7<J»)h  • n = (v  + 40 - (V  + 40h  , (5) 

where  i is  the  unit  positive  vector  along  the  x axis  (i  represents  the  unit  oncoming 
uniform  stream  equivalent  to  the  ship  speed  in  the  present  dimensionless  steady  flow 
formulation,  and  i + V4>  is  the  total  velocity  field)  as  defined  above,  N and  n are  the 
unit  inward  normal  vectors  to  (H)  and  (h),  respectively,  and  the  notation  ( )„  and 

Cl 

( )^  means  that  the  expression  between  the  parentheses  is  to  be  evaluated  on  (H)  or  (h) , 
respectively.  We  evidently  have  = 0 if  (h)  is  chosen  to  coincide  with  (H),  or  rather 
the  portion  of  (H)  below  the  plane  z = 0.  In  the  present  study,  it  is  convenient  to 
select  (h)  as  the  wetted  hull  of  the  ship  in  position  of  rest,  (H^)  say,  so  that  the 
fictitious  hull  (h)  = (R^)  and  the  real  hull  (H)  do  not  differ  much  (sinkage  and  trim 
are  not  very  significant  at  low  Froude  number),  and  q^  may  be  expected  to  be  small. 

It  is  shown  in  [1]  that  the  foregoing  "generalized  Neumann-Kelvin  problem"  can  be 
formulated  in  "integral  form"  given  by  the  integral  equation  (2.21),  which  here  becomes 


T 


i>*  + (<}>  - <j)*)G  da  + F2  4 [G(0<f>  + 

* Jh  n Tc  s 

r ? i 

Gq^dxdy  + Gq^da  - F ID  Gq^vyds 
'f  -'h  '’c 


[G(o<j>  + T<b  ) - (<J)  - $i.)G  ]yds  + 


where  the  significance  of  the  previously-undefined  symbols  will  now  be  explained.  The 
symbol  <p*  is  meant  for  <j>(x*)  where  x^  is  an  arbitrary  point  in  the  solution  domain  (d) 
including  its  boundary  (h)  + (f)  + (c)  , while  <J>  is  meant  for  $(x)  where  x represents 
the  "point  of  integration"  (integration  variable)  in  the  above  integrals;  the  point 
x*  thus  is  the  "field  point"  where  the  potential  is  being  evaluated,  while  x repre- 
sents the  "dummy"  variable  of  integration.  The  potential  <}>*  = $T(x*)  is  the  "initial 
potential"  given  by 


Vx*} 


G(x. ,x)v(x)da(x)  - F' 


G(x*,x 


)v  (s)y(s)ds 


The  function  G = G(x*,x)  in  equations  (6)  and  (7)  is  the  fundamental  solution 
(Green  function)  appropriate  for  the  problem,  that  is,  the  function  G(x*,x)  represents 
the  linearized  velocity  potential  for  the  disturbance  flow  caused  at  point  ^(x^.y^.z^O) 
by  a unit  "outflow"  at  point  x(x,y,z<0),  associated  with  a submerged  source  if  z < 0 
or  a flux  across  the  free  surface  if  z = 0,  in  an  oncoming  uniform  stream  along  the 
positive  x axis.  For  the  purpose  of  the  present  study  a convenient  expression  for 
the  Green  function  G is 

1 

4itG(x^,x)  “-7-pr  + ^'f  Im[exp(5)E1(?)  - ^]dt  + 

» F -1 

+ H(x')  ^ exp(— j T2)  sin(^r  I + tx)dt  , (8) 

F F F F 

—00 

• -►  . -►  -►  ^ ^ A 
where  r = |x*  - x|  is  the  distance  between  the  points  x^  and  x,  x = (x  ,y  ,z  ) = 

(x^  - x,  y*  - y,  z*  + z)  is  the  vector  joining  the  mirror  image  of  the  point  x with 

respect  to  the  plane  z = 0 to  the  "field  point"  x.,  r'  = |x'|  , £ = F 2(z^»/l  - t2  + 


y^t  + i|x"j)/l  - t , E^(£)  is  the  exponential  integral  defined  here  as  in  Abramowitz 
and  Stegun  [7],  H(x^)  is  the  Heaviside  step  function  defined  as  H(x^)  = 0 for  x"  < 0, 
that  is  upstream  from  the  point  x,  and  H(x')  = 1 for  x'  > 0,  and  T = J\  + t2 
Expression  (8)  for  the  Green  function  G may  readily  be  obtained  from  equations  (1), 

4» 

(2),  and  (6)'  in  Noblesse  [8].  The  above  expression  for  the  fundamental  solution 


+ 2 
Here,  we  used  the  relation  F 7r/r^=  - 


Im(l/£)dt,  which  may  easily  be  verified. 
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G will  actually  be  used  in  the  following  section  for  the  purpose  of  investigating  the 
"low-Froude-number  limit"  of  the  integral  equation  (6) . 

In  equations  (6)  and  '7) , the  symbol  da  in  the  surface  integrals  over  (h)  repre- 
sents the  differential  element  of  area  of  (h) , while  ds  in  the  line  integrals  around 
the  "waterline"  (c)  represents  the  differential  element  of  arc  length  cf  (c) . In 
th?.  line  integrals  around  (c) , we  have  V = V(s)  = n(s)  «i  where  n(s)  is  the  unit  inward 
normal  vector  to  (h)  at  point  s of  (c) , while  y is  defined  as  y = y(s)  = n'(s) *i,  where 
n"  is  the  unit  inward  normal  vector  to  (c)  in  the  plane  z = 0.  In  the  (fairly  common) 
case  when  the  surface  (h)  intersects  the  plane  z = 0 orthogonally,  we  have  n(s)  = n^(s) 
and  v(s)  = y(s) . In  the  first  line  integral  in  equation  (6),  the  symbols  a and  t are 
defined  as  a = s»i  and  T = t*i,  where  s = s(s)  is  the  unit  tangent  vector,  at  point  s, 
to  the  "waterline"  (c)  oriented  in  the  counterclockwise  direction  in  the  (x,y)  plane, 
and  t = t(s)  is  the  unit  tangent  vector  to  (h) , mutually  orthogonal  to  s(s)  and  the 
normal  n(s)  to  (h)  at  point  s of  (c) , and  pointing  downwards.  In  this  line  integral, 
the  notation  <j>  = 3<j>(s,  t,n) /8s , cj)  = 3<J)(s,t,n) /9t , and  G = 9G(x*,x)/9x  w.i>q  used  for 

shortness.  The  usual  notation  Gn  = VG(x*,x) »n(x)  was  also  used  in  the  first  surface 
integral  in  equation  (6).  It  will  finally  be  noted  that  the  axes  x,y,  and  z,  the 
"fictitious  hull  surface"  (h) , the  undisturbed  free  surface  (f),  the  "waterline"  (c) , 
the  elements  of  area  da  and  of  arc  length  ds,  and  the  unit  vectors,  i,  n,  n , s,  and 
t,  are  shown  in  figure  1. 

A method  of  solution  of  the  integral  equation  (6) , based  on  an  iterative  procedure 
starting  with  the  initial  approximation  ^(x^.)  given  by  formula  (7),  is  discussed  in 
some  detail  in  [1]  and  [2].  In  particular,  the  recurrence  relation  defining  the  kth 


approximation  <p 


in  the  iterative  sequence  of  approximations 


(1)  _ . (2)  .(3) 

i = 4>I , <p  , <p 


...  is  given  by  equation  (30b)  in  [2]  . The  formulas  defining  the  iterative  approxi- 
mations R,  to  the  wave  resistance  R associated  with  the  above-mentioned  approximations 
(k)  K 

<p  to  the  disturbance  velocity  potential  (p  may  also  be  found  in  [2].  The  main  object 
of  the  present  study  is  to  present  an  alternative  iterative  scheme  for  solving  the 
integral  equation  (6)  based  upon  the  assumption  that  the  value  of  the  Froude  number 
is  fairly  small,  which  is  in  fact  the  case  for  a large  class  of  commercial  ships.  Before 

Lf 

we  proceed  to  derive  these  alternative  low-Froude-number  iterative  approximations  <j), 

£f  k 

to  4>  , and  the  corresponding  approximations  R^  to  R,  we  begin  by  investigating  the 
limit  of  the  disturbance  velocity  potential  <J)  as  the  Froude  number  vanishes. 


In  fact,  the  recurrence  relation  (30b)  in  [2]  is  valid  for  k > 0,  as  it  is  indeed  in- 
dicated in  [2],  so  that  we  may  actually  regard  the  zeroth  approximation  = 0 as 

(k) 

the  initial  approximation  in  the  iterative  sequence  of  approximations  <j>  , k > 0 , 

which  thus  is  defined  by  4/°^  = 0,  = $ , <j/2\  ... 

I 
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2 . The  zero-Froude-number  limit 

In  the  differential  formulation  of  the  problem  of  steady  motion  of  a ship  stated  by 

1/2 

equations  (1)  through  (5)  the  Froude  number  F = U/(gL)  — where  L here  is  the  length 
of  the  ship  — appears  explicitly  only  in  the  free  surface  boundary  condition  defined  by 
equations  (2)  and  (4) . This  free-surface  condition  may  be  expressed  in  the  form 


% 

♦«  * - if 


(9) 


where  q^  is  defined  as 


’f  - k * f2{**k  + <IV*|2)X  + ^N|2>]  2 , ^ 1.^,2,  „ 

-lz  = - F (<)>x  + -^-l  VcJ)  I ) z = 0 


(9a) 


'V 


Equation  (9a)  shows  that  we  have  q^  -*■  0 as  F -*■  0,  so  that  the  free-surface  boundary 
condition  simply  becomes  the  "rigid-wall  condition"  <J>  = 0 in  the  "zero-Froude-number 

limit"  F = 0.  Physically,  as  the  Froude  number  vanishes,  the  force  due  to  gravity 
becomes  so  much  larger  than  inertia  forces  that  deformations  of  the  free  surface  are 
greatly  inhibited,  and  the  free  surface  indeed  appears  as  a rigid  wall,  at  least  on 
the  length  scale  of  the  ship. 

In  the  zero-Froude-number  limit  the  disturbance  velocity  potential  (p  becomes  the 
"zero-Froude-number  potential",  (p®  say,  which  thus  verifies  the  Laplace  equation 


V2<J>°  = 0 in  (d) 


(10) 


the  "rigid-wall  (Neumann)  condition" 


cj)  =0  on  (f)  , 

z 


(11) 


and  the  "Neumann  boundary  condition" 


- V = - n • i on  (h)  , 


(12) 


*! 


where  (h)  here  is  taken  as  the  wetted  hull  of  the  ship  in  position  of  rest,  i.e.  (h)  = 
(Hq)  + . The  "zero-Froude-number  potential"  may  readily  be  seen  to  be  identical  to  the 
disturbance  velocity  potential  of  the  flow  past  the  "double  ship",  consisting  of  (Hn) 


Sinkaga  and  trim  vanish  as  F -*■  0,  so  that  we  have  q^  = 0. 
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and  its  mirror  image  with  respect  to  the  plane  z = 0,  in  an  unbounded  fluid  (i.e.  with- 
out free  surface),  and  the  potential  <J>^  is  indeed  often  referred  to  as  the  "double- 
ship potential". 

It  is  also  interesting  to  examine  the  zero-Froude-number  limit  of  the  integral 

equation  (6),  for  this  provides  an  alternative  proof  of  the  fact  that  <p  -►  as  F -*■  0. 

For  this  purpose,  we  begin  by  examining  the  zero-Froude-number  limit  of  the  fundamental 

solution  (Green  function)  G(xA,x)  associated  with  the  Kelvin  free-surface  condition 
2 

G^  + F G^  = 0 on  z = 0,  and  given  by  formula  (8).  It  can  be  shown  [by  using  the 
asymptotic  expansion  of  the  exponential  integral  E^(£)  for  |?|»1]  that  the  first  in- 
tegral in  equation  (8),  which  represents  a nonoscillatory  near-field  disturbance  (as  it 

4 

is  shown  in  [8]),  is  0(F  ) as  F + 0,  while  the  last  integral,  which  is  associated  with 

ct  -l/F^ 

free-surface  gravity  waves,  is  0(F  e ) as  F 0,  with  a = 1 , 2 , or  2/3  depending  on 
whether  x^  is  inside,  outside,  or  exactly  on  the  boundary  of  the  famous  "Kelvin  V-shaped 
sector"  trailing  downstream  from  the  point  outflow  x.  We  then  have 

2 

G(x*,x)  = G°(x*,x)  + 0(F2)  + H(x')0(F6e~1/F  ) as  F + 0+  , (13) 


where  - 4/3  £B  = a-  2<_0,  and  G (x*,x)  is  the  "zero-Froude-number  Green  function" 
defined  by 


4ttG°(x.,x)  = - 1/r  - 1/r'  = - (1  + r/r')/r 


(14) 


The  function  G (x*,x)  defined  by  equation  (14)  clearly  is  the  Green  function  associated 
with  the  Neumann  condition  G^  = 0 on  z = 0,  and  thus  is  the  Green  function  appropriate 
for  the  zero-Froude-number  problem  stated  by  equations  (10)  through  (12).  More  pre- 
cisely, it  may  be  verified  that  the  fundamental  solution  G^(x.,x)  verifies  the  equations 


,2  0 


V G = <$(x  - x*)6(y  - y*)6(z  - z^)  in  z < 0 


G =0  on  z * 0 
z 


if  z*  < 0 , 


(15a) 

(15b) 


f 2 0 

'V  G =0  in  z < 0 


G ■ - <5(x  - x*)6(y  - y*)  on  z = 0 


if  z*  = 0 


(16a) 

(16b) 


where  6(  ) is  the  usual  Dirac  "5  function". 


It  ought  however  to  be  noted  that  r'  is  assumed  to  be  0(1)  as  F -*•  0 in  the  asymptotic 
approximation  (13),  which  is  therefore  not  valid  for  r'  = 0(F2);  in  other_,words , this 
"low-Froude -number  approximation"  is  not  uniformly  valid  with  respect  to  x'. 
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By  substituting  equation  (13)  into  equations  (6)  and  (7) , and  by  using  the  facts 


that  the  "iree-surf ace  flux"  q^  and  the  "hull  flux"  q^  are  0(F  ) as  F 


0,  it  may  be 

seen  that,  in  the  "low-Froude-number  limit",  the  integral  equation  (6)  becomes 


G Vda  + 


(H*)G  da  + 0(F  ) + 0(e 


-1/F 


) as  F -*■  0 


2 -1/F 

where  the  0(F  ) term  is  associated  with  a "near-field  disturbance"  while  the  0(e  ) 

term  corresponds  to  the  free-surface  gravity  waves,  which  thus  are  "exponentially  small" 

at  low-Froude-number . The  zero-Froude-number  limit  <j)^  of  the  disturbance  velocity 

potential  <p  thus  is  the  solution  of  the  "zero-Froude-number  integral  equation" 


[<f>°(x)  - <f>°(x*)  ]G°(xA,x)da(x) 


(17) 


where  the  "zero-Froude-number  initial  potential"  4>^.  Cx*)  is  defined  as 


+?<**>  = 


G (xJb,x)v(x)da(x) 


(13) 


It  may  be  verified  that  the  integral  equation  (17) , which  was  obtained  above  as  the 
zero-Froude-number  limit  of  the  integral  equation  (6),  is  in  fact  identical  to  the 
integral  equation  that  can  be  directly  derived  from  equations  (10)  through  (12)  and 
equations  (14)  through  (16)  by  applying  a classical  Green  identity  to  the  zero-Froude- 
number  potential  and  Green  function  G^  (in  the  manner  shown  in  Part  2 of  [1]), 
thereby  establishing  the  consistency  of  the  zero-Froude-number  limits  of  the  differential 
and  integral  formulations  of  the  problem  of  steady  motion  of  a ship  stated  by  equations 
(1)  through  (5)  and  equations  (6)  through  (8) , respectively. 

A straightforward  method  of  solution  of  the  integral  equation  (17)  consists  in 
using  an  iterative  procedure  based  on  the  recurrence  relation 


.0  /-*■  . ,(),-»■. 

<|)k+l(x*)  + 


r . 0 0 . , 0 /*^  ~^k 

[<f>k(x)  - 4>k(x*)  ]Gn(x*,x)da(x) 


k > 0 


(19) 


which  immediately  follows  from  equation  (17) , and  beginning  with  the  "zeroth  approxi- 


mation" <J>q(x*)  = 0. 


This  "zeroth  approximation"  <}>q  becomes  exact  in  the  limiting  case 


of  a "knife-blade-like"  or  "needle-like"  ship-hull  form,  hich  suggests  that  the  above 
iterative  procedure  corresponds  to  a "thin-,  flat-,  or  slender-ship  perturbation  method 
of  solution".  In  particular,  the  "first  approximation"  in  the  sequence  of  iterative 
approximations  <J>k,  k >_  0,  defined  by  the  recurrence  relation  (19)  and  the  "zeroth  ap- 
proximation" <t»Q  - 0 is  given  by  = <J>®,  so  that  the  "initial  potential"  defined 


*~3 
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by  formula  (18)  may  thus  be  regarded  as  a first-order  zero-Froude-number  slender- 
or  thin-  , or  flat-  ship  approximation. 

As  a matter  of  fact,  the  initial  potential  in  some  respects  appears  to  gener- 
alize the  classical  first-order  thin-  and  slender-body  perturbation  approximations,  as 
it  will  now  be  shown.  For  convenience,  we  begin  by  expressing  the  initial  potential 
given  by  formula  (13)  in  the  equivalent  form 


i ° * -1 

<|,I(x*)  =4? 


n(x)  • i da(x) 

h"  [ (x*-x) 2 + (y*-y)2  + (z*-z)2]1/2 


where  equation  (14)  was  used,  and  (h")  represents  the  "double-hull  surface"  consisting 
of  the  surface  (h)  and  its  mirror  image  with  respect  to  the  plane  z = 0.  With  the 
coordinates  x and  x^  made  dimensionless  in  terms  of  a reference  length  taken  as  the 
length  of  the  ship,  the  dimensionless  length  of  the  double  hull  (h')  is  of  course  equal 
to  unity,  so  that  in  the  case  of  a thin  ship  form  we  have  jy|«l,  and  expression  (20) 
may  be  approximated  as 


. -1  f n(x)  • i da(x) 

' 4?  L “ *2""  2 “ , TTTJT 

Jh  L(x.-x)  + y.  + (z.-z)  ] 


for  |y^.|  ■ 0(1),  that  is  in  the  "outer  (thin-ship)  limit"  in  the  usual  language  of  the 
"method  of  matched  asymptotic  expansions"  (see  for  instance.  Van  Dyke  [ 9 ]).  If  the 
equation  of  (h')  is  expressed  by  means  of  the  two  equations  y = b+(x,z)  and  y = - b (x,z) 
corresponding  to  the  two  sides  of  the  surface  OO  , the  surface  integral  (21)  may  be 
transformed  into  the  following  double  integral  over  the  projection,  (lO  say,  of  (h") 
onto  the  plane  y = 0 


4>j(x*)  * 


*L 


t^(x,z)dxdz 

x2  ~ 2 , , *2,1/2 

[ (x*-x)  + y*  + (z*_z)  J 


where  t(x,z)  = b (x,z)  + b (x,z)  is  the  "thickness"  of  the  hull.  In  the  usual  case  of 
a ship  with  port  and  starboard  symmetry,  we  have  b~(x,z)  = b+(x,z),  so  that  t(x,z)  = 
2b(x,z),  and  expression  (22)  becomes  identical  to  the  classical  first-order  perturbation 
approximation  for  a thin  body.  This  well-known  thin-ship  (body)  perturbation  approxi- 
mation may  thus  be  seen  to  correspond  to  the  thin-ship  limit  of  the  initial  potential 
defined  by  formula  (18),  in  the  same  way,  of  course,  that  the  classical  Mrchell 
thin-ship  approximation  corresponds  to  the  thin-ship  limit  of  the  initial  potential 
given  by  formula  (7),  as  it  was  shown  in  [1],  Indeed,  the  recurrence  relation  (19) 
proposed  for  solving  the  zero-Froude-number  integral  equation  (17)  is  essentially 


.J 
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identical  to  the  recurrence  relation  (2.28a)  in  [1]  related  to  the  integral  equation  (6) 
furthermore,  both  of  these  recurrence  relations  are  associated  with  the  use  of  the 
approximations  = 0 and  $ = 0 as  "zeroth  approximations",  which  clearly  corresponds 
to  a perturbation  method  of  solution  appropriate  for  thin-  or  slender-ship  forms,  as 
it  was  noted  previously. 

In  the  (usual)  case  of  ship  forms  which  are  not  only  thin  but  also  flat,  i.e. 
slender,  we  have  |z|«l  as  well  as  |y|«l,  and  expression  (22)  for  the  initial  potential 
may  be  further  approximated  as 


•fa  ' £ 


\(x)dx 

r , ^2  2 2.1/2 
t (x*-x)  + y*  + z*] 


(23) 


where  A(x)  is  given  by 


X(x) 


• d (x) 

t (x,z)dz 
~d (x)  X 


(23a) 


with  d(x)  representing  the  local  (at  section  x)  draft  of  the  ship  hull  (h) . For  usual 

ship  forms  we  have  X(x)  = da(x)/dx,  where  a(x)  is  the  cross-sectional  area  of  the  double 

ship  form  (h'’).  The  approximation  (23),  which  was  obtained  above  as  the  "slender-ship 

2 2 1/2  2 2 1/2 

outer  limit"  [based  on  the  assumption  (y  + z ) « 1 with  (yy.  + z^)  = 0(1)]  of 

expression  (20)  for  the  initial  potential  4>^,  may  thus  be  seen  to  be  identical  to  the 
"outer  solution"  in  the  well-known  first-order  slender-ship  approximation  based  on  the 
use  of  the  method  of  matched  asymptotic  expansions. 

It  may  also  be  interesting  to  examine  the  "slender-ship  (body)  inner  limit"  of  the 
surface  integral  (20)  defining  the  initial  potential  <J>  . For  this  purpose,  it  is 
convenient  to  formally  introduce  the  "inner  variables"  n,  C and  n*,  C*  related  to  the 
"outer  variables"  y,  z and  yA,  z*  by  means  of  the  equations  y = en,  z * ££,  and  y*  “ er.fc 
z^  * where  e is  the  slenderness  parameter  attached  to  the  ship  (e  may  be  taken  as 

the  ratio  of  the  largest  transversal  dimension  of  the  ship  to  its  length).  In  terms  of 
these  inner  variables,  formula  (20)  may  be  expressed  in  the  form 


4)l(x**n*»^*) 


■>  -¥■  ■+■  •+  9 —1/9 

n(x,eri,eg)  »i  [1  - (n*i)  ] eda 

[(x*-x)2  + e2((n*-n)2  + (c*-02>]1/2 


where  (a)  is  the  intersection  curve  of  the  double-hull  surface  (h')  with  any  plane 
x » constant,  and  eda  is  the  differential  element  of  arc  length  of  (a).  By  performing 
the  change  of  variable  x * x*  - e£,  we  may  obtain 
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e 


nCx^-eC.en.^S)  *T  [l-(n*i)2]_1//2 

[C2  + (n*-n)2  + U*-02]1/2 


cda 


In  the  limit  e -*■  0 — and  for  ea  <_  1 - £a  with  0 < a < 1,  that  is  away  from  the 

ship  ends  (bow  and  stern)  x^  a 0 and  x^  * 1 — we  may  then  obtain 


0 .,1  n(xA,en,eC)  *i  f°° 

»!<*., I,.?.)  ' ">»  j1_(J.£)2 jl/2 


d£_ 


[C2  + (n*-n)2  + (C*-02]1/2 


where  the  longitudinal  slope  n»i  clearly  was  assumed  to  be  continuous.  The  inner 
integral  in  the  above  expression  is  divergent;  however,  proper  interpretation  of  this 
integral  yields 


« 1 _ f / n ^ nCx^.en ,££) *i 

in/(n*-n)z  + (S*-cr  '[1~^~I)2]i/2  edG  * (24) 

2 2 1/2 

This  then  shows  that  in  the  "inner  region"  (y^  + z*)  - 0(e)  — and  away  from  the 

ship  ends  x*  * 0 and  x*  * 1 — the  initial  potential  <J>  is  locally  two  dimensional, 

0 1 

with  x*  appearing  as  a parameter,  i.e.  we  have  (n*,C*;x^) , which  is  a well-known, 

and  indeed  intuitively  self-evident  property  of  flows  about  slender  bodies.  More 

precisely,  formula  (24)  expresses  the  inner  approximation  to  the  initial  potential 

$j(n*,C*;xA)  as  the  potential  due  to  a distribution  of  two-dimensional  sources,  of 

^ ^ 2 l/2 

strength  (n«i) / [l-(n*i)  ] , along  the  intersection  curve  (a)  of  the  double-hull 

surface  (h")  with  the  plane  x = x. . This  inner  approximation  to  the  initial  potential 
thus  is  clearly  not  identical  to  the  inner  solution  in  the  classical  slender-body 
perturbation  approximation,  which  is  also  quasi  two-dimensional  of  course,  but 
satisfies  the  condition  of  no  flow  across  the  "framelines"  (a)  exactly,  while  this 
condition  is  only  approximately  verified  by  the  initial  potential  <J>^.  In  any  case,  we 
must  evidently  keep  in  mind  that  — in  the  present  "longitudinal"  motion  — the  most 
significant  part  of  the  ship  hull  surface  in  fact  are  its  ends  (bow  and  stern) , where 
flow  deceleration  and  acceleration  mainly  take  place,  rather  than  the  approximately 
uniform  "middle-body"  section  to  which  the  above  inner  approximation  is  restricted. 

It  may  be  interesting  to  incidentally  indicate  here  an  alternative  form  of  the 
recurrence  relation  (19)  . This  alternative  form  can  be  obtained  by  noting  that  the 
kth  approximation  <f>^  on  the  right  side  of  equation  (19)  verifies  the  following  identity 


G°(x*,x)<J>£n(x)da(x)  + [4>°(x)  - <t>°(xy.)  ]G°(xA,x)da(x)  , 

' h ' n 
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which  can  be  established  by  applying  a classical  Green  identity  to  the  potential  <t*?(x) 

0 f k 

and  the  Green  function  G (x^,x)  . Use  of  the  above  identity  and  of  equation  (18)  in 

equation  (19)  then  yields  the  following  alternative  recurrence  relation 

* + [ G°(**>*Hv(x)  (x)]da(x>  | , k ^ 0 , (25) 

I ih I 

which  must  also  be  associated  with  the  zeroth  approximation  $q(xa)  = 0.  While  the 
recurrence  relation  (19)  is  actually  better  suited  than  the  (equivalent)  alternative 
recurrence  relation  (25)  for  practical  purposes  of  numerical  calculations,  the  latter 
relation  readily  lends  itself  to  an  appealing  physical  interpretation,  namely  one  keeps 
distributing  sources  of  strength  V + = (i  + V<*>^) *n  equal  to  the  fluid  flux  across  (h) 

kn  o k 0 

associated  with  each  successive  approximation  <j>^  until  the  zero-flux  (v  + = 0)  hull 

boundary  condition  (12)  is  verified  within  the  desired  accuracy,  which  then  provides  a 
physical  insight  into  the  iterative  scheme  associated  with  the  recurrence  relation  (19) ; 
in  short,  this  iterative  scheme  thus  might  be  interpreted  as  a successive  "leakage 
stopping"  scheme. 

It  is  instructive  to  consider  the  simple  particular  case  when  the  ship  hull  surface 
(h)  is  the  "bottom  half"  of  a triaxial  ellipsoid  with  principal  axes  oriented  parallel 
to  the  x,  y,  and  z axes  and  corresponding  principal  dimensions  2a,  2b,  and  2c,  respec- 
tively, so  that  the  length,  beam,  and  draft  of  this  simple  ship  hull  form  are  2a,  2b, 
and  c,  respectively.  As  it  was  noted  previously,  the  velocity  potential  of  the  dis- 
turbance flow  due  to  a hull  form  (h)  in  the  zero-Froude-number  limit  is  identical  to 
the  potential  of  the  flow  due  to  the  "double  hull"  (h')  consisting  of  (h)  and  its  mirror 
image  with  respect  to  the  plane  z = 0 in  an  unbounded  fluid  (i.e  without  a free  surface). 
In  the  particular  case  of  a triaxial  ellipsoid  now  being  considered,  an  exact  expression 
for  the  disturbance  velocity  potential  is  in  fact  known;  this  expression  (see  for 
instance  Havelock  [10])  is 


'(x*)  - 


-1 


2rr(2  - a)  3x4 


dv(x) 


lx*  - xl 


(26) 


where  V-  represents  the  interior  volume  of  the  ellipsoid,  i.e.  V-  is  the  volume  bounded 
externally  by  the  surface  (h')  of  the  ellipsoid,  dv  is  the  differential  element  of 
volume,  and  a is  a constant  defined  by  the  integral 


This  expression  may  readily  be  verified  by  comparing  it  with  the  integral  equation 
(17),  in  which  one  needs  only  replace  <p  by  and,  in  formula  (18),  V by  (-  4^)  as 
it  may  be  seen  from  equation  (12) . 
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a * abc 


r 

;0 


(A  + a2)"3/2[(A  + b2)(X  + c2)]'1/2dX. 


(26a) 


t * . . — 


Formula  (26)  may  be  expressed  in  the  form 


\ -1  f 3 , 1 

* (**>  ‘ 2tt( 2 - a)  J_  3x<  _ 


-)dv(x) 


X,.  - X 


2ttC2  - a) 


±(- 

3xV 


-)dv(x)  . 


The  latter  volume  integral  may  finally  be  transformed  into  the  surface  integral 


j.0 ,-*■  \ -1  f n(x)  »i  , 

* <x*)  Jh-  l^Tif daU)  ' 


(27) 


where  n(x)  is  the  unit  inward  normal  vector  to  (h")  at  point  x,  in  accordance  with  the 
definition  introduced  previously. 

Comparison  of  expressions  (27)  and  (20)  then  shows  that  the  exact  potential  cj5*^  in 
fact  is  proportional  to  the  initial  potential  <p®,  that  is  we  have 


(28) 


$0(x  ) * k I 

L _ 1 r_  _) 


where  the  constant  of  proportionality  k is  given  by  k 3 1/(1  - a/2).  By  performing  the 

2 

change  of  variable  A * a t in  the  integral  (26a) , and  upon  introducing  the  notation 
8 = b/a  (B  = beam/length  ratio)  and  6 = c/2a  (6  = draft/length  ratio),  we  may  obtain 

I k(M)  - l/fl  - 86  | (t  + 1)~3/2  [(t  + B2)  (t  + 462)]'1/2  dtl  I . (28a) 

! !\ o 

It  may  be  seen  that  k(8,<$)  -+■  1 if  B -*■  0 (thin-ship  limit),  or  if  6 0 (flat-ship  limit) 

or  if  both  8 0 and  6 -►  0 (slender-ship  limit),  in  agreement  with  the  fact  that  the 

initial  potential  must  become  exact  in  the  thin-,  flat-,  or  slender-ship  limits. 
Curves  k(B;<5)  are  represented  in  figure  2 for  0 B 1 and  0 5 1/2  (the  ellipsoid 

becomes  a sphere  for  8*1,  5 ■ 1/2).  It  may  be  seen  that  for  8 “ -2  and  6 * .1,  which 
may  be  regarded  as  typical  values  of  the  beam/length  and  draft/length  ratios,  we  have 
k = 1.06,  so  that  the  initial  potential  4>^  given  by  formula  (18)  in  fact  provides  a 
fairly  accurate  approximation  to  the  velocity  potential  of  the  zero-Froude-number 
disturbance  flow  due  to  an  ellipsoidal  hull  form  of  typical  ship-like  dimensions. 

To  be  sure,  real  ship-hull  forms  are  not  ellipsoidal  in  shape.  Nonetheless,  the 
initial  potential  seems  likely  to  provide  a realistic  approximation  to  the  zero- 
Froude-number  disturbance  potential  , and  at  least  ought  to  provide  a fairly  good 
first  approximation  for  solving  the  integral  equation  (17)  iteratively  on  the  basis 
of  the  recurrence  relation  (19).  Furthermore,  the  remarkable  result  expressed  by 


— — ,,  --^i- 1 jj— arr  - 1 -||j 


J 
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equation  (28)  suggests  that  the  "modified  initial  potential"  given  by  k^,  where  k(8,6) 
is  the  function  defined  by  formula  (28a)  and  depicted  in  figure  2 — with  8 and  6 taken 
as  the  beam/length  and  draft/length  ratios,  respectively  — may  provide  a slight  improve- 
ment of  the  original  initial  potential  <J>j  given  by  formula  (18)  . 

Equation  (28)  also  suggests  that  a promising  alternative  approach  to  solving  the 
integral  equation  (17)  consists  in  expressing  the  solution  <f>^  of  this  integral  equation 
in  the  form  4>^(x*)  ■ k(x*)<f>^(xA)  . By  substituting  this  relation  into  equation  (17), 
and  expressing  the  term  <J>^(x)  - <{>^(x*)  = k(x)<fc!j!(x)  - k(x^)<{> j(x^)  in  the  form 
k(xA)[d>j(x)  - ^(x^)]  + [k(x)  - k(xA)  ](pj(x)  , we  may  obtain 


k(x<()<|i°(xjk)  * + k(x*) 

L 1 Jh  L 


(x)  - $>j(x*)  ]G°(x*,x)da(x)  + 


+ 1 [k(x)  - k(x*)  ]4>I(x)Gn(xJlt,x)da(x) 

' h 

By  multiplying  both  members  of  this  equation  by  4>^ ( xA)  , using  the  relations  <f>^(x*)  = 
k(x*)<j>!j!(x*)  and  4^(x)  * k(x)<J>!j*(x) , and  rearranging,  we  may  finally  obtain  the  following 
alternative  form  of  the  integral  equation  (17) 


I 4>  (x*)  « + 


f (<t>°(x)<t>j(xA)  - 4>°(x*)<j>j(x)  ]G^(xJ|(,x)da(x) 

‘ h 

4»°(x*)  - [4>j(x)  - <}>!j!(x*)  ]G°(x*,x)da(x) 

' h 


where  the  potential  $j.(x^),  which  corresponds  to  the  initial  potential  associated  with 
the  modified  integral  equation  (29),  is  given  by 


I ^<x*)  ■ [♦£(**)  )2/{4>5<x*)  - | [4>j(x)  - 4>?(V]G°(x*,x)da(x)|  I 


(29a) 


The  above  approximate  expression  for  the  zero-Froude-number  potential  (x*)  is  exact 
in  the  case  of  a triaxial  ellipsoid  in  translation  parallel  to  its  major  axis,  as  it 
was  shown  previously,  and  thus  seems  likely  to  provide  a fairly  good  approximation  in 
the  case  of  real  ship-hull  forms.  It  may  be  interesting  to  note  that  by  performing  a 
formal  binomial  expansion  (and  retaining  only  the  first  two  terms  in  this  expansion) 
of  expression  (29a)  based  on  the  assumption  that  the  surface  integral  in  the  denominator 
is  small  compared  with  the  potential  4>j(x*),  we  may  obtain 

$°(x*)  * $°(x*)  + (4>j(x)  - <^(x*)  ]G°(x*,x)da(x)  , 


J 


*“0  ^ 0 v 0 

which  shows  that  <j>T(x*)  - 4u(x*)>  where  (J)9(xjJc)  is  the  second  approximation  in  the  se- 

quence  of  iterative  approximations  4>^(x^) , k 0,  defined  by  the  recurrence  relation  (19) 

It  may  incidentally  be  noted  that  the  above  approach  can  obviously  also  be  applied 

to  the  Neumann-Kelvin  integral  equation  (6).  Thus,  by  expressing  the  solution  <{>(x*)  of 

the  "linearized"  Neumann-Kelvin  integral  equation  (6)  — in  which  the  nonlinear  free- 

surface  flux  q^  and  the  hull  flux  q^  (associated  with  effects  of  sinkage  and  trim)  are 

ignored  — in  the  form  4>(x*)  = k(xA)<frj(xA) , where  ^(x*)  is  the  initial  potential  given 

by  formula  (7),  we  may  obtain  the  following  expression  for  the  approximation  iL.(x*) 

-0  1 

corresponding  to  expression  (29a)  for  the  approximation  <p^(x*) 


L_. 


V**)  = 


{♦I  - j 

V Jy 


($X  - ^)Gnda  ~ 


*1 


[G(<x4>J  + T4»J)  - (<P1  - <fri)G  Jyds 

S L w X 


where  the  notation  <£*  = 0^. (x^)  and  <p^  = <j>^(x)  was  used  for  shortness.  This  approxi- 
mation will  be  examined  further  in  Part  5 of  this  study. 

It  may  finally  be  noted  that  while  the  integral  equation  (17)  and  the  modified 
integral  equation  (29)  are  valid  for  x*  on  the  hull  surface  (h)  and  in  the  flow  domain 
(d) , and  in  principle  may  thus  be  used  for  determining  in  (h)  + (d) , it  may  be  more 
expedient  in  practice  to  use  the  integral  equations  (17)  or  (29)  for  determining  <t>^  on 
(h),  and  then  to  obtain  <p^  in  (d)» strictly  outside  (h) , by  means  of  the  well-known 
expression 


(30) 
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3.  First-order  low-Froude-number  wave-resistance  formulas 

As  it  was  explained  in  Part  3 of  this  study  [2],  the  term  <p  - <}>*  in  the  integrands 
of  the  first  and  second  integrals  on  the  right  side  of  the  integral  equation  (6)  [which 
is  valid  for  x*  on  the  hull  surface  (h)  + (c)  and  in  the  flow  domain  (d)  + (f) ] may  be 
replaced  by  (j)  in  the  far  field  of  the  ship,  that  is  for  x*  -*>  <*>  , since  we  have  -►  0 
as  |xj  -*•<»,  and  therefore  |<{>  J « |<j)|  for  x*  in  the  far  field  and  x in  the  near  field. 
Furthermore  , insofar  as  one  is  only  interested  in  the  surface  waves  trailing  far  down- 
stream from  the  ship,  the  first  three  terms  in  expression  (8)  for  the  Green  function 
G(x*,x),  which  represent  a nonoscillatory  near-field  disturbance,  may  be  neglected. 

The  disturbance  velocity  potential  <[>(x*)  far  downstream  from  the  ship  may  then  be  ap- 
proximated as 


4>(x  ) 'u  F2Re  ^ 

* IT 


u+tV/2 


exp  j-|(l+t2)  + if-4  + e)  X+tHdt 


which  expresses  the  potential  in  terms  of  a familiar  superposition  of  elementary 
plane  waves. 

The  function  Q(t)  in  the  above  expression  is  the  "Kochin  free-wave  spectrum 
function"^  defined  as 


' 

C2(t)  = ftg(t)  + 4>Enda  + [E(a<J)g  + T4»t)  - <)>Ex]yds  + 

'h  ' c 


+ Eq^dxdy  + Eq^da  “ F -r  Eq^vyds  , ( 
' f ■'h  ' c 

where  the  function  ftg(t),  which  is  referred  to  as  the  "zeroth  approximation  to  the 
Kochin  function"  £2(t)  since  it  corresponds  to  taking  <J>  = 0 as  approximation  to  the 
disturbance  potential  <p  in  the  ship  near  field,  is  given  by 


Vfc> 


E(x;t)v(x)da(x)  - F‘ 


t E(‘; 


t)v  (s)y(s)ds 


with  the  function  E(x;t)  defined  as 


E(x;t) 


(l+t2)3/2 


exp{'%(1+t2)  * iiJS2+ 

Vp  v p ) 


rIt  will  be  noted  that  there  is  here  a slight  difference  in  the  definition  of  the  Kochin 

function  used  in  the  present  Part  4 and  in  Part  3 [2]  of  this  study;  specifically,  the 

Kochin  function  f2(t)  defined  in  this  Part  4 is  related  to  the  Kochin  function  ft(9)  of 

3 

Part  3 by  the  equation  f2(t)  = Q(9)sec  9. 


The  dimensionless  wave  resistance,  R say,  defined  as  R = R / pU  L where  R is  dimen 
sional  (p  is  the  density  of  water) , can  then  be  determined  from  the  above-defined 
Kochin  spectrum  function  £)(t)  by  means  of  the  well-known  "Havelock  wave-resistance 
formula" 


|fl(t) |2  (1  + t2)"5/2  dt 


Further  details  on  the  derivation  of  the  foregoing  formulas  may  be  found  in  Part  3 of 
this  study  [2], 

The  discussion  of  the  previous  section  immediately  suggests  that  a "first-order 

•dF 

low-Froude-number  approximation",  (t)  say,  to  the  Kochin  function  ft(t)  may  be  ob- 
tained by  approximating  the  disturbance  potential  $ in  the  ship  near  field,  in  formula 
(31),  by  the  potential,  \\?  say,  which  may  be  taken  as  either  the  "exact  zero-Froude- 
number  potential"  <J>^  given  by  the  solution  of  the  "zero-Froude-number  integral  equa- 
tion" (17)  or  the  "zero-Froude-number  slender-ship  approximation"  k(B,6)<J>*j!  given  by 
formulas  (18)  and  (28a)  , [The  approximation  defined  by  formula  (29a)  will  not  be 


used  here.]  The  first-low-Froude-number  approximation  (t)  then  is  given  by 


— ■ — — 

jximation  fifF(t) 


^ (t)  = nQ(t)  + 


ip°Enda  + F2  [E(cnj/^  + x^P_)  - -^E^yds  + 


Eq^dxdy  + Eq^da  ~ F2  Eq^\ 

* f Jr* 


where  the  hull  flux  q^  is  defined  by  formula  (5),  with  <p  replaced  by  that  is  we  have 


%• (v  + - (u  + *°)h  ’ 


(36a) 


while  the  nonlinear  free-surface  flux  q is  defined  by  formula  (4) , with  <|>  also  re- 
0 * 

placed  by  ^ of  course.  By  performing  a Taylor  series  expansion  of  the  right  side  of 
equation  (4)  about  the  plane  z = 0,  we  may  obtain 

qf  = F2  j(|V<i>|2)x  + |?<l>-V|V<f>|2  - (<Px  + y|V4>|2)<f>2z}  + 0(F4)  , 

where  the  expression  on  the  right  side  is  to  be  evaluated  at  the  plane  z = 0.  Replace- 
ment of  <p  by  , and  use  of  equations  (38a,b)  below,  then  yield 


“ F2  + +i 


°|2  + |V2^°|2V2^°  1 + 0(F4)  , 


(36b) 


where  is  the  two-dimensional  differential  operator  (8^,3^).  By  using  formula  (32) 


for  the  zeroth  approximation  (L(t)  in  expression  (35),  we  may  express  the  first  low- 


tiF , 


Froude-number  approximation  (2^  (t)  in  the  form 


«fF(t) 


[E(V  + q^)  + ip^E  ]da  + 
Jh  h n 


Eq^dxdy 


n 


i 

- F2*k  [E{v(v  + q°)  - (cn|;°  + t^°)}  + 4)°Ex]uds  I . 

j 


(37 


If, 


An  interesting  alternative  expression  for  the  above-defined  approximation  (2.  (t) 

0 -*■  1 

may  be  obtained  by  applying  a classical  Green  identity  to  the  functions  \p  (x)  and 


E(x;t)  in  the  domain  (d) , as  it  will  now  be  shown.  We  begin  by  noticing  that  the 


.0 


zero-Froude-number  potential  ip  (x)  verifies  the  following  equations 


vV  = 0 in  (d)  , 


(38 


0 on  (f)  , 


= - V + (v  + (j/*)  on  (h)  , 

n n 


j^|  't-  1/  | x | ^ as  |x|  -*■  00  , 


(38 

(38 

(38 


where  we  actually  have  V + ip^  = 0 if  ip®  is  taken  as  the  exact  zero-Froude-number 
0 n 

potential  $ , as  it  may  be  seen  from  equation  (12) . We  also  note  that  the  function 


E(x;t)  defined  by  formula  (33)  verifies  the  equations 


V^E  =0  in  z < 0 , 


E + F E = 0 in  z < 0 , 

Z XX  — 


-4  2 2/ 2 2 

E ^ F (1  + t)  “ exp[z(l  + t )/F  ] as  z -*■  - 


(39 

(39 

(39 


as  it  may  readily  be  verified. 

It  follows  from  equations  (38a, d)  and  equations  (39a, c)  that  the  following 
Green  identity 


°E  - E^°jd 
n n^ 


a 
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holds.  Equations  (38b,c)  and  equation  (39b)  then  yield 


t p^E  da  + f Evda  - f E(v  + ip^)da  - F2  ip^E  dxdy  = 0 

h n Jh  Jh  n if  xx 


By  expressing  the  term  iJ/^E^  in  the  last  integral  in  the  orm  lj ®Evv  - (lP^^x  ~ Bf>x)x  + 


Eip  , and  by  using  the  well-known  Green  identity 
XX 


fy°Ex  - &l£)x  dxdy  = -^  (i^°Ex  - Ey°)dy  = C^°Ex  - E’J^Uds  , 


where  we  used  the  relation  dy  = - yds  along  (c)  [with  y = n'*i  as  it  was  defined  pre- 
viously in  connection  with  equations  (6)  and  (7)],  we  can  obtain 


Evda  + ip 
h ^h 


'E^da  + F2  tD  (E^x  - ij/^Jyds  = E(v  + da  + F 
' c •'h 


2 ri0 

Eip  dxdy  . 
J xx 


As  it  was  shown  in  Parts  2 and  3 of  this  study  [1,2]  we  have  = cnf^  + + v^  , 

0002  0 xstn 

which  becomes  \p^  = cnp^  + - V + v(v  + ijj^)  by  virtue  of  equation  (38c)  . We  may 

then  obtain 


Evda  - F2  Ev2yds  + j ip^E^da  + F2  tp  [ E (oip ® + Tip®)  - i p^E  ]yds  = 
c J h n J c S C x 

F2  E\p0  dxdy  + E (v  + da  - F 2 E (v  + i[j^)vyds 

h xx  Jh  n h n 


By  substituting  the  above  equation  into  expression  (35) , and  using  formulas  (32)  and 
(36a),  we  can  finally  obtain  the  following  alternative  expression  for  the  first  low- 

£f 

Froude-number  approximation  (t) 

] fi^F(t)  = F2  E(^x  + ^f)dxdy  + E(v  + Rda  " f2^  E(v  + OHVyds  ' ’ (40) 

1 _Jh | 

where  q^  is  defined  as  q^  = q^/F2. 

With  (h)  taken  as  the  wetted  hull  of  the  ship  in  position  of  rest,  and  ip^  chosen 
as  the  exact  zero-Froude-number  potential  < , we  have  (v  + ij/j)  = 0 by  virtue  of  equa- 
tion (12) , and  expression  (AO)  becomes 

0^F(t)  =■  F2  E^^  + 9f]dxdy  . (41) 

If,  furthermore,  the  nonlinear  free-surface  flux  q^  is  neglected,  we  obtain  the 
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Guevel-Kayo  low-Froude-number  approximation",  ft  (t)  say,  given  by 

i ~r  : i 

| ft~  (t)  = F E(x,y,z=0;t)4>  (x,y,z=0)dxdy  | ; (42) 

|_G^ Jj “ 

indeed,  the  above  approximation  to  the  Kochin  function  ft(t)  is  essentially  the  same  as 
the  low-Froude-number  approximation  which  was  first  obtained  by  Guevel,  Vaussy  and 
Kobus  [ 3 ] , and  is  precisely  identical  to  the  formula  later  derived  by  Kayo  [ 4 ] . As 
it  may  be  seen  from  equation  (37),  an  alternative  expression  for  the  Guevel-Kayo 


approximation  is 


^GK(t) 


(Ev  + 4>GE  )da  - E{v2  - (cr4>G  + x<j>G 

n Tc  L st 


+ *°Ex} 


(42a) 


which  may  in  fact  be  equally  convenient  as  formula  (40)  for  purposes  of  numerical 
calculations . 

Approximation  (41) , in  which  the  nonlinear  free-surface  flux  qG  is  retained,  may 
be  shown  to  be  identical  to  the  "Baba-Maruo  low-Froude-number  approximation",  ft^(t) 


say,  which  thus  is  given  by 


= + Eq^dxdy  = F2  e(4>!L  + 5^)  dxdy  | » 


f J Jf  -^xx 


(43a, b) 


where  qG  = F2qG  is  defined  by  formula  (36b),  with  ip G replaced  by  . Indeed,  it  can 
£ t _0 

easily  be  verified  that  the  term  + q^  may  be  expressed  in  the  form 


■4 + ■ [o + 0 + tiv0'2)]  + [*X + Hv°i2)_  _ 


(43c) 


which  is  identical  to  expression  (11)  in  Baba  [ 5 ] and  expression  (32)  in  Maruo  [ 6 ] . 

If  the  potential  <p®  is  taken  as  the  approximation  given  by  k(|3,<5)4>G,  where  k(3,<5) 
is  the  function  of  the  beam/length  ratio  8 and  draft/length  ratio  6 defined  by  formula 

(28a)  and  <{)  is  the  initial  potential  given  by  equation  (18),  an  alternative  form  of 

1 ZF 

expressions  (37)  and  (40)  for  the  first-order  low-Froude-number  approximation  ft  (t) 

0 £ 

may  be  obtained  by  exploiting  the  fact  that  the  potential  ip  now  is  defined  not  only 
in  the  flow  domain  (d)  but  also  in  the  "interior  domain"  (d^) , i.e.  in  the  domain 
bounded  externally  by  the  hull  surface  (h)  and  the  portion,  (f^)  say,  of  the  plane 


z = 0 inside  the  "waterline"  (c) . In  fact,  it  may  easily  be  shown  that  the  potential 
ip®  verifies  the  equations 
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V2/  = 0 in  (d±)  , 


(44a) 


t[i  =0  on  (f . ) , 

z i 


(44b) 


ipO  = kv  + k e<p®  on  (h)  , 

rn  In 


(44c) 


where  the  fact  that  the  potential  ^ = kcf;^  is  associated  with  a distribution  of  sources 
of  density  kv  on  the  surface  (h)  was  used  in  equation  (44c)  [specifically,  we  have 

- eip®  = kv,  where  the  superscripts  i and  e refer  to  the  "interior"  and  "exterior" 

n n e 0 
sides  of  (h) , respectively],  and  the  superscript  e in  the  symbol  <j>  is  meant  to 

0 in 

clearly  indicate  that  the  normal  derivative  <J>T  of  the  zero-Froude -number  initial 
0 

potential  ij>  is  to  be  evaluated  on  the  exterior  side  of  the  hull  surface  (h)  . 

By  virtue  of  equations  (39a)  and  (44a),  the  following  Green  identity 


fy°En  - &l£)da  = (<P°E„  ~ Eip°)dxdy 


Z Z' 


holds.  By  using  equations  (44b),  (39b),  and  (44c),  we  may  obtain 

f 

\p®E  da  - k Evda  - k E e<J>^  da  + F^  ip^E  dxdy  = 0 
i n , , In  ~ xx 

h f . 

1 

By  expressing  the  term  ip^E  in  the  form  (ip^E  - Etp^)  + E^  , and  using  the  Green 

XX  X XX  XX 

identity 

0°Ex  - &l£)  dxdy  = t£>  Cip°Ex  - E^x)  dy  = - 1£»  fy°Ex  - yds  , 

J f ) r* 


where  we  used  the  relation  dy  = - yds  along  (c) , we  can  obtain 

ij Pe  da  - k Evda  - k E e<j)^  da  + F^  *4  (E<p^  - ip®E  )yds  + F^  Eip^  dxdv  = 0 

h n 'h  Jh  ln  h x x h.  xx 

1 

As  it  was  noted  previously,  we  have  , which  becomes  = 

0 0 2,  e 0 ^ s c n x 

cnp^  + + kv  + k by  virtue  of  equation  (44c)  , so  that  we  may  obtain 

tp®E  da  + F^  [Efcnj/*  + Tt|^)  - t|^E  ]pds  = kS)n(t)  - F^  Eip®  dxdy  + 

, n I st  x U _ xx 

'h  ' c ' f ^ 

+ k E e<J)jnda  - F^k  t£>  E e<J>^avyds  , 

^h  ' c 


where  formula  (32)  was  used.  By  substituting  the  above  equation  into  expression  (35) 

ZT 

for  the  first-order  low-Froude-number  approximation  ft^  (t),  and  using  equations  (36a) 
and  (32).  we  may  then  obtain  the  following  alternative  expression  for  the  first-order 
low-Froude-number  slender-ship  approximation,  f^Fs^(t)  , say 


-iFs l 


(t)  = kftQ(t)  + 


E (v  + 4;°]  H da  - F2  41  E (v  + Hvyds 

h ' n 


- F 


f . 

l 


xx 


Ety  dxdy  + Eq  dxdy 


With  the  potential  ^ explicitly  expressed  as  k(3,<5)4>^  we  finally  have 


^Fs£(t)  = kftQ(t)  + j 


E(v  + k°()>^)Hda  - F 2 E(v  + k°4>F) Hvyds 


- 


L 


r 0 . 1 , . , _2.  2 
E (f^dxdy  + F k 


r 0 1 , , 
t q^dxdy 


J 


(45) 


where  the  zero-Froude-number  initial  potential  <J>  , given  by  formula  (18),  was  written 
in  the  form  for  convenience,  and  the  nonlinear  free-surface  flux  ^qF  is  defined  as 


0 I 


•V2|V2^|2  + I Vl!2v2^ 


(45a) 


as  it  may  be  obtained  from  equation  (36b) . The  first  two  terms  on  the  right  side  of 
equation  (45a)  are  0 ( | \7<{>!j!  | 2)  , while  the  terms  between  brackets  are  o(|V<}>^|^j,  so  that 
these  "third-order"  terms  might  be  neglected  for  usual  slender  ship-hull  forms. 

The  above  "first  low-Froude-number  approximation"  fil^FsL  may  be  regarded  as  an 
approximation,  specifically  a low-Froude-number  approximation,  to  the  "first  slender- 
ship  approximation"  ft.  defined  by  formula  (22)  in  Part  3 of  this  study  [2].  Indeed, 
comparison  of  formula  (45)  above  and  of  formula  (22)  in  Part  3 shows  that,  besides  the 
fact  that  the  multiplicative  correction  function  k(8,<5)  was  not  introduced  in  expression 
(22)  in  Part  3 (i.e.  we  simply  have  k = 1 in  this  formula),  the  present  low-Froude- 
number  approximation  (45)  essentially  corresponds  to  approximating  the  initial 
potential  <J>^,  given  by  formula  (7),  by  the  zero-Froude-number  initial  potential  <p® 
given  by  formula  (18)  . In  addition,  however,  the  above  expression  (45)  contains  a 
surface  integral  over  the  ship  "waterplane"  (f^)  that  has  no  counterpart  in  expression 
(22)  in  Part  3.  The  numerical  simplification  resulting  from  the  substitution  of  the 
potential  <p ^ given  by  formula  (7)  by  the  zero-Froude-number  potential  <p®  given  by 
formula  (18),  which  mainly  stems  from  the  much-simpler  form  of  the  zero-Froude-number 


Green  function  G given  by  formula  (14)  in  comparison  with  the  Kelvin  Green  function 

G given  by  formula  (8),  is  significant,  and  indeed  a main  recommendation  of  the  low- 

Froude-number  slender-ship  approximation  (45)  in  comparison  with  the  slender-ship 

approximation  (22)  of  Part  3 resides  in  the  fact  that  the  approximation  (45)  does  not 

require  evaluation  of  the  fairly  complicated  Kelvin  Green  function  (14) . The  approxi- 
£Fs-£ 

mation  ft^  defined  by  formula  (45)  also  appears  to  be  somewhat  simpler  than  the  low- 
Froude-number  approximations  ft^  and  ft^  defined  by  formulas  (42, a)  and  (43a, b)  in 
that  these  approximations,  based  on  the  exact  zero-Froude-number  potential  evidently 
require  the  preliminary  determination  of  the  potential  <j)^,  that  is  the  preliminary 
solution  of  a zero-Froude-number  integral  equation  (the  method  of  Hess  and  Smith  [11] 
is  usually  used  for  this  purpose) . 

/Fs  £ 

Expression  (45)  for  the  first  low-Froude-number  slender-ship  approximation  ft^  (t) 
may  be  written  in  the  form 


(t)  = nf!s£(t) 


E q.dxdy 


where  the  linearized  first  low-Froude-number  slender-ship  approximation  ft^  (t)  is 
given  by 


ftf*s£(t)  = EAda  - F2t£>  EAvyds  - F2k  E dxdy  i , 

'h  Tc  XX 

i I 


with  the  source  density  A defined  as 

A = kv  + (v  + k °<}>*)H  , ( 

n n 

as  it  can  be  obtained  by  using  formula  (32) . An  interesting  alternative  expression 
for  the  approximation  is 


(Ev  + k^E  ) da  - F2  [E{v2  - kfa0^1  + + k (p^E  ]yds 

hx  n / s t l x 

J c 


which  can  readily  be  obtained  by  substituting  k(f)  for  4)  into  expression  (42a)  . The 

lYsi  1 

above-defined  approximation  ft^£  is  the  simplest  of  the  low-Froude-number 

approximations  obtained  in  this  section.  As  a matter  of  fact,  comparison  of  formula 

(32)  for  the  zeroth  approximation  ft  with  formulas  (48)  and  (47)  for  the  approximation 
£Fs£  U 

ft£^  shows  that,  from  the  computational  point  of  view,  the  latter  approximation 
requires  a relatively  simple  modification  of  the  basic  zeroth  approximation  ftA(t). 


In  summary,  several  (evidently  related)  low-Froude-number  approximations  to  the 
Kochin  free-wave  spectrum  function  ft(t)  have  been  obtained.  The  common  characteristic 
feature  of  these  various  first-order  low-Froude-number  approximations  resides  in  that 
they  are  based  upon  approximating  the  disturbance  velocity  potential  (j>  in  the  ship 
near  field,  in  formula  (31),  by  the  "zero-Froude-number  potential"  t which  was 
taken  as  either  the  "exact  zero-Froude  number  potential"  <J>^,  given  by  the  solution  of 
the  "zero-Froude-number  integral  equation"  (17),  or  as  the  "zero-Froude-number  slender- 
ship  approximation"  k(0,<$)4>T  defined  by  formulas  (18)  and  (28a).  While  the  first-order 
low-Froude-number  approximations  ^ and  are  somewhat  more  "sophisticated"  than 

the  basic  "zeroth  approximation"  given  by  formula  (32) , which  is  simply  associated 
with  the  approximation  t)>  = 0 in  the  ship  near  field,  these  approximations  admittedly 
are  still  fairly  "crude".  More  "refined"  (second-order  low-Froude-number)  approxi- 
mations  (^  and  ) can  however  be  defined  by  using  a more  "realistic"  approxi- 

mation to  the  potential  <J>  in  the  ship  near  field,  in  formula  (31),  than  the  "wave- 
free"  zero-Froude-number  potential  ip  . The  analysis  of  section  2 readily  suggests 

that  natural  approximations  for  <J>  are  the  "first-order  low-Froude-number  approximations" 
£f  5/Fs& 

<£,  and  4>  say,  defined  by  evaluating  the  various  unknown  terms  on  the  right  side  of 

0 

the  integral  equation  (6)  on  the  basis  of  the  zero-Froude-number  potential  \p  (taken  as 
<\P  or  as  k<$Kj!)  as  an  approximation  to  <p  . These  first  low-Froude-number  approximations 
to  the  potential  <J>  will  now  be  determined. 
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First-order  low-Froude-number  approximations  to  the  velocity  potential 


As  it  was  noted  above,  the  first-order  low-Froude  number  approximation  to  the 
potential  <f>  may  be  defined  by  approximating  $ in  the  unknown  terms  on  the  right  side  of 
the  integral  equation  (6)  by  the  zero-Froude-number  potential  ip^ . We  thus  have 

$f(x*)  » ^(x*)  + 0|»°  - ^)Gnda  + F24  [G(cnp°  + - (<J/°  - 4°)Gx]uds  + 

'h  -*c 

+ f Gq^  dxdy  + Gq^  da  - F2  Gq^  vyds  , (49) 


where  the  potential  4>T(x..)  is  the  initial  potential  defined  by  formula  (7),  and  q, 

0 i n 

and  q^  are  the  hull  flux  and  the  nonlinear  free-surface  flux  given  by  formulas  (36a) 

and  (36b),  respectively.  By  using  equation  (7),  expression  (49)  for  the  approximation 

tF 

may  be  written  in  the  form 

| $fF(x.v)  = [ [G(v  + q®)  + (tp°  - ^)G  ]da  + Gq°f  dxdy 

j 1 Jh  h if  f | 

| - F2  j>  [G{v(v  + q°)  - (cnp°  + ti|»°)}  + 0J/°  - *J)G  ]uds  ! . (50) 

! 1C J 

An  interesting  alternative  expression  for  the  first-order  low-Froude-number  ap- 

£F  0 

proximation  <(>^  can  be  obtained  by  exploiting  the  fact  that  the  potential  ip  verifies 

the  equations 


V^tp®  =0  in  (d)  , 


,0  . „2 .0  _2 . 0 

4 + F 4 = F 4 on  (f)  , 

z xx  xx  ’ 


^ - [V  - (v  + ^)h]  on  (h)  , 


(51a) 


(51b) 


(51c) 


as  it  can  readily  be  seen  from  equations  (38a, b,c),  so  that  ip  must  satisfy  the 

4. 

following  identity  1 

4'*  =■  G[V  - (v  + \p°),  ]da  - F2*£  G[v  - (v  + ip0)  ]vyds  - F2  ( G^°  dxdy  + 

Jh  nh  Jc  nh  Jf  xx 

+ f (<P°  - 'jfrc  da  + F2i>  tG(Olp°  + xifi?)  - (4°  - 4>J)G  ]Uds 

Jh  11  Jr  S L X 


(51d) 


Equation  (51d)  may  be  written  down  at  once  by  comparing  equations  (51b, c),  (2),  and 

(3),  which  show  that  one  simply  needs  to  replace  the  terms  , qf,  and  (v  + q ) by  the 

0 2 0 0 * h 
terms  , - F 4^.  and  [v  - (v  + 4^)^],  respectively,  in  the  integral  equation  (6). 
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where  we  used  the  notation 


Ax.)  and  i Ax). 


and  (36a)  into  equation  (49),  we  may  obtain  the  following  alternative  expression  for 

£f 

the  approximation  <J>^ 


| <x*> 


1 


Ax.)  + F2 


G(C+^f)dxdy  + 


L_. 


G(v+AHda  - F2*D  G(vV)  vudsl, 
n n I n n 

__C I 


where  qG  = qG/F2  as  it  was  defined  previously. 


If  the  fictitious  hull  (h)  is  taken  as  the  wetted  hull  of  the  ship  in  position  of 


rest  [so  that  we  have  q^  = 0 in  expression  (50),  by  virtue  of  equation  (36a)],  the 
potential  is  taken  as  the  exact  zero-Froude-number  potential  <}>G  [so  that,  by  virtue 


of  equation  (12),  we  have  (v  + ip  ) = 0 in  expression  (52)],  and  the  nonlinear  free- 

0 _ 2-0  n “ 

surface  flux  q^  = F q^  is  neglected,  the  potential  <J>^  becomes  the  "Guevel-Kayo 


potential",  <{>  say,  given  by  the  alternative  expressions 

(jK 


L_ 


A**)  + F2 


I 


G(x.,x,y,0)<f>  (x,y,0)dxdy 


xx 


(53) 


I 


*GK(X*) 


.0  x0 


[Gv+(4>  -<j)*)G  ]da  - F‘ 


* 
J c 


"1 


[G{v2-(a<J>G+r<])G) } + (0°-4>°)Gx]yds 


J 


corresponding  to  expressions  (52)  and  (50),  respectively.  If  the  nonlinear  free- 

0 ZF 

surface  flux  q^  is  retained,  the  potential  i}^  becomes  the  "Baba-Maruo  potential". 


^BM  Say’  given  bF 


^BM^X*^ 


*GK(x*>  + 


L_ 


G(^xx  + ^f)dxdF 


(53b, c) 


„0 


If  the  zero-Froude-number  potential  ^ is  taken  as  the  slender-ship  approximation 
k(B,<5)<t>G,  where  k(B,6)  is  the  function  of  the  beam/length  ratio  3 and  draft/length 
ratio  6 defined  by  equation  (28a)  and  is  the  potential  given  by  formula  (18),  an 


alternative  form  of  expressions  (50)  and  (52)  for  the  first  low-Froude-number  poten- 

£F  0 0 

tial  <j>^  can  be  obtained  by  exploiting  the  fact  that  the  potential  = kcj) ^ now  is 


defined  not  only  in  the  flow  domain  (d)  but  also  in  the  interior  domain  (d.)  bounded 

1 0 

by  the  hull  surface  (h)  and  the  "waterplane"  (f^).  More  precisely,  the  potential  ip 


verifies  the  following  equations 


n2 . 0 
V ip 


in  (di) 


By  using  equations  (51d) , (7), 


(52) 


(53a) 


(54a) 


J 
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1 


<i>°  + fV  = F2k  °4>I  on  (f  ) , 

z xx  xx  i 


ty®  * kv  + k et()^  on  (h)  , 
n In 


(54b) 

(54c) 


as  it  may  be  obtained  from  equations  (44a, b,c).  It  can  then  be  shown  that  ip  must 
satisfy  the  following  identity 


(4>°  - n£)Gnda  + F24>  [G(oi|£  + Ti|>°)  - OP0  - 1[|J)G  luds  = 
J c 


G(kv  + k e<J>?  ) da  - F2  4>  G(kv  + k e<j>?  )vyds  - F2k 
h In  h In 


f . 
1 


G ^4>^xdxdy 


(54d) 


t 


By  substituting  equation  (54d)  into  equation  (49),  and  using  formulas  (7)  and 


(36a) , we  may  then  obtain  the  following  expression  for  the  "first-order  low-Froude- 

£Fs£ 

number  slender-ship  approximation",  <J>  say. 


n 


Fs£ 


♦r  (x*)  = fefrjU*) 


r 


+ J, 


0,1 


G(v  + k 9*)Rda  - Ffc  *p  G(v  + k V^)nvyds 


t 


0,1, 


n 


rn'H 


2 1 0 1 2 2 1 0 1 
- F k G u4>  dxdy  + F4k  G dxdy 

J f xx  J f f 

i 


(55) 


where  the  zero-Froude-number  initial  potential  <(>  , given  by  formula  (18),  was  written 
in  the  form  ^<J>^  for  convenience,  and  the  nonlinear  free-surface  flux  is  given  by 


formula  (45a)  . Equation  (55)  may  be  written  in  the  form 


,1 Fs£, 

*1 


,£Fs£, 


2,  2 


♦,(**>  = (**}  + F k 


G °q*  dxdy  | , 


(56) 


tYs'l 


where  the  linearized  low-Froude-number  slender-ship  potential  ^ ' can  [by  using 


formula  (7)]  be  expressed  as 


,£Fs  l 


hi  (X*} 


GXda  - F‘ 


G> 

* n 


Xvyds  - F k 


~7  i 

G <J)  dxdv  | , 
xx  y I 


(57) 


i > j 


with  the  source  density  X given  by  formula  (47a)  . An  interesting  alternative  expression 


iFsZ 

for  the  approximation  is 


The  identity  (54d)  can  be  established  by  applying  the  approach  used  in  Part  2 of  this 
study  [1],  with  straightforward  modifications  accounting  for  the  fact  that  equations 
(54a, b,c)  correspond  to  an  interior  (rather  than  exterior)  problem. 
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,£Fs£,-*-  . 

hi  (x*> 


[G\H-k(V-V)G  ]da  - F2i  [ G { v2 -k ( a°(J> ) }+k(V-V) G ]yds  ,(58) 

* n l st  x x 


0 1 0 

which  can  be  obtained  by  substituting  k <j>  = k<J>j  for 


into  expression  (53a) . 


Comparison  of  expression  (55)  for  the  first-order  low-Froude-number  slender-ship 
. .IFsl  . . , „ , , . . , 


approximation  d> . 


and  expression  (2.27)  in  Part  2 of  this  study  [1]  for  the  second- 


£Fs  £ 

order  slender-ship  approximation  i}^  shows  that  the  potential  can  be  regarded  as 

a low-Froude-number  approximation  to  the  potential  Specifically,  it  may  be  seen 

that,  besides  the  fact  that  the  multiplicative  correction  function  k(S,6)  was  not  intro- 

JlFsJl 

duced  in  expression  (2.27)  for  the  approximation  ^ , the  potential  essentially 

corresponds  to  approximating  the  initial  potential  $ given  by  formula  (7)  by  the 

0 ^ 

zero-Froude-number  initial  potential  <f>  given  by  formula  (18),  although  an  additional 

1 iFsl 

integral,  over  the  ship  "waterplane"  (f  ) , appears  in  expression  (55)  for  <}>  . A 

1 iTsl  i 

main  recommendation  of  the  low-Froude-number  approximation  , in  comparison  with 

the  approximation  obviously  resides  in  the  notable  numerical  simplification 

associated  with  the  substitution  of  the  potential  <(>!?  for  the  potential  4>  . As  a matter 

^ iFsl  ^ 

of  fact,  from  the  computational  point  of  view,  the  approximation  4>^  , expressed  in  the 
form  given  by  formula  (58)  say,  corresponds  to  a relatively  simple  modification  of  the 
basic  initial  potential  $ , given  by  formula  (7),  associated  with  the  approximation 

4>  = 0 in  the  unknown  terms  on  the  right  side  of  the  integral  equation  (6).  Actually, 

JlFsJl 

the  first  low-Froude-number  slender-ship  approximation  4>  , associated  with  the 

0 1 

approximation  4>  * k<f>^  in  the  unknown  terms  in  the  integral  equation  (6)  , may  just  as 
well  be  regarded  as  a low-Froude-number  approximation  of  the  second  slender-ship  ap- 
proximation <J>2.  as  it  was  just  discussed,  or  as  a generalization  of  the  initial 
potential  (fi^  ( = first-order  slender-ship  approximation  j>^)  , and  it  may  indeed  readily 
be  verified  that  expressions  (55)  through  (58)  become  expression  (7)  if  k is  replaced 
by  zero.  In  conclusion  to  this  section,  it  will  be  noted  that  expressions  (50),  (52), 
(53,a,b,c),  (55),  (56),  (57)  and  (58)  for  the  several  (obviously  related)  first  low- 
Froude-number  approximations  to  the  potential  <{>  clearly  correspond  to  expressions  (37), 
(40),  (42, a),  (43a, b),  (45),  (46),  (47),  and  (48),  respectively,  for  the  first  low- 
Froude-number  approximations  to  the  Kochin  free-wave  spectrum  function  0,  and  indeed 
the  latter  approximations  to  the  Kochin  spectrum  function  could  evidently  have  been 
derived  from  the  approximations  for  the  potential  <p,  that  is  following  the  present 
section  rather  than  preceding  it. 
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5 . Second-order  low-Froude-number  approximations 


tY 

A second-order  low-Froude-number  approximation,  say,  to  the  disturbance 

potential  <p  may  readily  be  defined  by  approximating  the  potential  p on  the  right  side 

of  the  integral  equation  (6)  by  the  previously-determined  first  low-Froude-number 

XY 


approximation  <f> 


We  may  then  obtain 


XY 


*2  (X*}  = Vx*>  + 


ZY  ZY  2 1 ZY  ZY  ZY  ZY 

(^1  -^l*)Gnda  + FZ"b  [G(C<|>~4t4>")  - «>" -4>i*)G  luds  + 
h ' c 


ZY  Zy  2 X Zy 

Gq~dxdy  + Gqjida  - Y**p  Gq^vyds  , 

f Vi  ' r* 


(59) 


£f  Zy 

where  the  nonlinear  free-surface  flux  q ..  and  the  hull  flux  q are  given  bv  formulas 

ti  £F  hi 

(4)  and  (5),  respectively,  in  which  <p  is  replaced  bv  . By  using  formula  (7)  for 

,-*■  ZY  £ 

the  initial  potential  <p^(x*)  , expression  (59)  for  (x*)  may  be  written  in  the  form 

ZY 


XY.+  \ 
<Py  (x*) 


h + v>  + <«?  - C,GnIda  + \{ 


Jf  G^fidxd^ 

F2^  [G{V(V  + q^)  - <0$**  + T^)}  + (<^F  - <0^)Gx]ud£ 


(60) 


An  alternative  expression  for  the  above  second-order  low-Froude-number  approxi- 

Zy  Zy 

mation  <J) ^ can  be  obtained  by  exploiting  the  fact  that  the  first  approximation  $ 
verifies  the  equations 


2 ZY 

7 ^ «0  in  (d) 


XY  V2XY 
<P,  + F 1 

lz  lxx 


pln 


- qf  on  (f) 


Xy. 


[v  - (v  + 4>ln)h]  on  (h) 


(61a) 

(61b) 

(61c) 


as  it  can  readily  be  seen,  so  that  we  have  the  identity 


X* 

h* 


j 


G[v  - (v^)h]da  ~ F2 fj?  G[v  - (v^F)h]vuds  + Gq^dxdy  + 
h * c ' f 


(4>f -^*)Gnda  + F2 
h • c 


lj>  [G(o^F+t^F)  - (t)>fF-())ff)G.Juds 


'1*'  xJ 


(61d) 


^F  i.F  tF  Z.F 

where  <J>^  = (x^)  and  (J»^n  = VfJ)^  *n.  By  using  equations  (61d) , (5),  and  (7)  into 
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equation  (59) , we  may  then  obtain  the  following  alternative  expression  for  the  approxi- 

£F 

mation 


J.F.-+  x JX.-+  , 

<t>2  (x*)  = <{)1  (x*)  + 


G(v'H>ln)Hda  - G(v+<J)^)Hvyds  + 

h ^ c 


■If 

Gtt^  dxdy 


(62) 


£f 


0 


£F  tF 

where  the  free-surface  flux  it  is  defined  as  tt 

1 ^F  ^ L tF 

(4),  in  which  <p  is  replaced  by  4>-^  * and  equation  (61b)  that  vr^  is  given  by 


qfl  " qf V 


It  follows  from  equation 


IF  rtF  v2(.lF  r . „ . £F  , 2-,  lv.lF  _ | v.tF , 2T| 

'l  1*1.  + F K*  + (IV*,  I )x  + 2Vh  •'='1^1  I )J2  . f2(^F  + I|70tt|2)  • 


(62a) 


From  the  alternative  expressions  (60)  and  (62)  for  the  second  low-Froude-number 
tF 

approximation  $ to  the  potential  <p  , we  can  easily  derive  the  following  alternative 

^ tF 

expressions  for  the  second-order  low-Froude-number  approximation  ft^  to  the  Kochin 

free-wave  spectrum  function  ft(t) 


tF 

(t) 


- F 


[E(V  + q^F)  + ^FEn]da  + 


tF 

Eqfldxdy 


f 

J c 


[E{v(v  + q^)  - (a<^F  + T<t>^)}  + 4>^FEx]yds 


Pv 

E tt^  dxdy  | 


(63a) 


(63b) 


Second-order  low-Froude-number  approximations  R_  to  the  wave  resistance  R can  then  be 

1 tF 

readily  defined  by  using  the  above  second-order  approximation  to  the  Kochin  free- 

wave  spectrum  function  ft  in  the  Havelock  wave  resistance  formula  (34). 

£Fs£,  tFsl  tFsH 

Second-order  low-Froude-number  slender-ship  approximations  <)>  , > anc* 

to  the  potential  $,  the  Kochin  spectrum  function  ft,  and  the  wave  resistance  R,  respect- 

tF  tFsZ 

ively,  can  similarly  be  defined  by  merely  replacing  by  in  expressions  (60), 

(62),  and  (63a, b) . Of  main  practical  interest  may  be  the  second-order  low-Froude- 

tFsl 

number  slender-ship  wave  resistance  approximation  R associated  with  the  approxi- 

IFsl  2 

mation  ft  to  the  Kochin  spectrum  function  ft,  corresponding  to  the  use  of  the 

^ £.Fs 

linearized  first-order  slender-ship  approximation  ' defined  by  formula  (58)  as  an 


approximation  to  the  potential  <p  in  formula  (31)  for  0. 
in  equation  (63a)  readily  yields 


A slight  change  of  notation 
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re,  -L.  £Fs ls  ^ ,£Fs lr  , , , r lYsl  . , 

„ [E(V  + qhU  } + En]da  + . EqfU  dxdy 

h r 


W2JL  rrr  . . £Fs l.  . J-Vsl.,  . ,l7slr:  , . „ 

- F ^ [E{v(v  + qhl£  ) - (a<j)l£>s  + T^^)}  + 4>l£  Ex]lids  , (e 

where  the  notation  4>  />sE  = 3<})^s^/3s,  and  = 3<|)^S^/3t  was  used,  the  hull  flux 

ZFsl  i'c'»s  Lt  i-t>t  lv  si 

is  given  by  equation  (5)  with  <j>  replaced  by  4>  ^ , that  is  we  have 


(v  + 4C1  - (V  + 


lZ,n'H 


l£,n'h  ’ 


(64a) 


and  the  nonlinear  free-surface  flux  q » is  defined  by  formula  (4)  in  which  $ must 

£Fs£  1 

evidently  also  be  replaced  by  <j>  ^ ; by  performing  a Taylor  series  expansion  of  the 

right  side  of  equation  (4)  about  the  plane  z = 0,  and  retaining  only  the  terms  that 
r2  2 

are  0 (f  | V(f>  | ) as  a first  approximation,  we  may  obtain 


iFsl  A „2  rnn.£Fs£|2,  .ZFsl  XFsll 
qfU  - F L(l  1 ' ♦u.x  4’u,zzj 


(64b) 


where  the  expression  on  the  right  side  is  to  be  evaluated  on  (f) , i.e.  on  the  plane 
z = 0 . 

It  is  interesting  to  compare  the  above  second-order  low-Froude-number  slender- 

iFsZ 

ship  approximation  and  the  straightforward  first-order  slender-ship  approximation 

given  by  equation  (18)  in  Part  3 of  this  study  [2].  In  terms  of  the  notation  used 
in  the  present  Part  4,  this  slender-ship  approximation  is  given  by  the  expression 


S3  (t)  = [E(v  + q*)  + < pJE  ]da  + Eq^dxdy 
i jh  n n ) £ 1 

- F2^)  [E{v(v  + qj)  - (o<{>^  + T<J>*)}  + <j)IEx]uds 


where  = <p_  is  the  initial  potential  defined  by  formula  (7) , the  hull  flux  qj'  is 
I h 

given  by  equation  (5)  with  <J>  replaced  by  4> ^ , that  is  we  have 


% ’ (v  + *n>H  - (v  + *X  • 


(65a) 


and  the  nonlinear  free-surface  flux  is  given  by  equation  (3a)  in  [2],  which  becomes 


<f  ■ p2  [d7»ii2)x  - + p2Oz]  • 


(65b) 
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Expressions  (64, a, b)  and  (65, a, b)  clearly  are  identical,  except  for  the  fact 

that  the  initial  potential  in  expressions  (65)  for  Q1  becomes  the  linearized 

1 IF si1 

first  low-Froude-number  slender-ship  approximation  4>_ » in  expressions  (64)  for 

iFsl  i/C 

(2  . Any  difference,  notably  from  the  computational  point  of  view,  between  the 

^ ZFsZ  ZFsZ 

approximations  £1^  and  ^ must  then  stem  from  the  potentials  <J>j-  and  <(i  ^ , which 

are  given  by  formulas  (7)  and  (58),  respectively.  As  it  was  already  noted  in  the 

,ZFsZ 


previous  section,  the  problem  of  evaluating  the  potential  $ 


1 Z 


numerically  is  not 


significantly  more  difficult  than  that  of  evaluating  the  potential  <p  , so  that  the 

ZFsZ  1 ZFsZ 

wave-resistance  approximations  and  associated  with  the  approximations 

and  (2^  are  comparable  from  the  point  of  view  of  required  numerical  calculations . 
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6 . Alternative  derivations  of  the  low-Froude-number  approximations 

Z.F  £Fs£  £F  ZFsZ 

The  low-Froude-number  approximations  <f)^  , <p^  and  <)> ^ , <p^  were  obtained  in 
the  previous  sections  as  first  and  second  approximations  in  a sequence  of  iterative 
approximations  associated  with  the  integral  equation  (6)  and  the  "zeroth  approximation" 
ipO , which  was  taken  as  the  "exact"  zero-Froude-number  potential  given  by  the  solution 
of  the  zero-Froude-number  integral  equation  (17)  or  as  the  zero-Froude-number  slender- 
ship  approximation  k. (S , 6) cf>^  defined  by  formulas  (18)  and  (28a).  These  low-Froude- 
number  approximations  can  also  be  obtained  by  using  two  related  alternative  approaches 
based  on  formally  expressing  the  potential  <5  in  the  form 


cj)  = ip®  + \p 


(66) 


either  in  equations  (1)  through  (5)  or  in  equation  (6) , that  is  in  the  differential  or 
integral  formulation  of  the  problem.  These  alternative  approaches  will  now  be  briefly 
examined . 

By  substituting  equation  (66)  into  equations  (1)  through  (5) , which  must  evidently 
be  satisfied  by  the  potential  <j)  = ip  + ip,  and  using  equations  (38a, b,c)  verified  by  the 
zero-Froude-number  potential  \\P , we  may  obtain  the  following  equations  to  be  satisfied 
by  the  "low-Froude-number  potential"  ^ 


V‘  ip  « 0 in  (d) 


(67a) 


K + F\*  - - (f2*L  + 4 + <;>  “*  <f>  • 


>„  ■ - [(x 


^ = - |(v  + t^n)H  + 0|>n)H  - (lp„)J  on  (h)  , 


Vh] 


(67b) 

(67c) 


where  the  nonlinear  free-surface  flux  qf  + q'  is  identical  to  q (that  is  we  have 

o ^ 0 

q^  + q^  = q^)  given  by  formula  (4)  in  which  4>  is  replaced  by  ip  +!(;  , with  q^.  defined  as 


q° = f z + p2{4 + n^°i2)x  + 

and  q'  defined  by  the  relation 


* - - F2(*°  +I|V*°|2) 


f2[*°J  • (67d> 

z=0 


qf  = qf  - qf  • 


(67e) 
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In  addition  to  equations  (67a, b,c),  the  usual  radiation  condition  of  no  waves  upstream 
from  the  ship  must  evidently  be  satisfied.  It  may  be  verified,  by  performing  a Taylor 
series  expansion  of  the  right  side  of  equation  (67d) , that  the  low-Froude-number  ap- 
proximation to  the  nonlinear  free-surface  flux  qG  is  given  by  equation  (36b)  obtained 
previously.  It  will  also  be  noted  that  the  expression  on  the  right  side  of  the  hull 
boundary  condition  (67c)  vanishes  if  the  fictitious  hull  (h)  is  taken  as  the  wetted 


hull  of  the  ship  in  position  of  rest  and  effects  of  sinkage  and  trim  are  neglected, 


and  if  is  taken  as  the  exact  zero-Froude-number  potential  , since  we  then  have 


(v  + 4>  ) = 0 by  virtue  of  equation  (12)  . 

n H. 


The  "generalized  Neumann-Kelvin  problem"  defined  by  equations  (67a, b,c)  and  the 
radiation  condition  may  be  stated  in  integral  form,  that  is  in  the  form  of  an  integral 
equation,  in  the  manner  shown  in  Part  2 of  this  study  [1].  Specifically,  the  following 
integral  equation  can  be  obtained  for  the  low-Froude-number  correction  potential  ip 


>k**)  = jh  {K)h  + (Oh]da  - ?2f  g[K)h  + yl]v,jds  + 

+ jf  G^2^L+qf+qf)dxdy  + j (^-^JGnda  + [G(^t)  - 0M>J<yyds 


(68) 


where  the  notation  fib  = fib  - fib  V was  used  for  shortness,  and  U;.  and  ip  are  meant 

n'h  v n-'H  v n- h * 


for  ip(x+)  and  \p(x) , respectively,  with  x*  representing  an  arbitrary  point  in  the 


domain  (d)  + (f)  + (h)  + (c)  while  x is  the  "integration  point",  in  accordance  with 
previous  practice.  The  integral  equation  (68)  may  be  expressed  in  the  form 


iKx.)  = iMxJ  + 


X v 


- ^*)Gnda  + F2  |g (c^g  + T^t)  - ('P  ~ 'J'jcjjuds 
+ | GfyjJda  - F2^  GfojJJvyds  + f Gq'dxdy  , 


+ 


(69) 


where  the  potential  ^(x*)  is  defined  as 


VX*}  = 


G(F2^  + q°)dxdy  + G(v  + Hda  - F2  4 g(v  + i^^vyds  . (69a) 


t, 


The  integral  equation  (68)  may  readily  be  obtained  by  comparing  equations  (67b, c)  and 


equations  (2)  and  (3) , which  show  that  one  simply  needs  to  replace  the  terms  q^  and 


(v  + 9^]  by  the  terms  + qG  + qf)  and  [ fv  + + (<l>  )j|]»  resPectivedF»  in  the 

integral  equation  (6) . 
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The  above  potential  ij;  (x^) , which  clearly  regroups  the  various  "known  terms"  (that 
is,  the  terms  involving  the  presumably-known  zero-Froude-number  potential  p alone)  in 
the  integral  equation  (68)  provides  a natural  initial  approximation  [corresponding  to 
the  neglect  of  the  various  "unknown  terms",  i.e.  the  integrals  involving  the  unknown 
potential  ip,  in  the  integral  equation  (68)]  for  solving  this  integral  equation  itera- 
tively, in  the  manner  used  previously  in  Part  2 of  this  study  [1].  The  potential  <p 

defined  by  substituting  the  initial  approximation  \p  given  by  formula  (69a)  into  ex- 

L ZF 

pression  (66)  may  be  seen  to  be  identical  to  the  approximation  <}>..  obtained  previously 

ZF^ 

in  equation  (52) . The  second  low-Froude-number  approximation  <p ^ could  likewise  be  re- 
derived here  as  the  second  approximation  in  the  sequence  of  iterative  approximations 
associated  with  the  integral  equation  (69) . 

The  integral  equation  (69)  for  the  low-Froude-number  correction  potential  was 
derived  above  from  the  "differential  formulation"  of  the  problem,  specifically  from 
equations  (1)  through  (5),  into  which  equation  (66)  was  substituted.  As  one  would 
expect,  the  integral  equation  (69)  can  also  be  obtained  by  substituting  equation  (66) 
directly  into  the  "integral  formulation"  of  the  problem,  that  is  the  integral  equation 
(6),  and  it  may  be  interesting  to  briefly  consider  this  alternative  approach  here.  By 
substituting  equation  (66)  into  the  integral  equation  (6),  expressing  the  nonlinear 
free-surface  flux  q^  and  the  hull  flux  q^  in  the  forms  q^  = q^  + q',  as  defined  by 
equations  (67d,e),  and  qh  = q°  + 6^)^  , with  q°  = (v  + - (u  + as  defined  by 

equation  (36a)  and  (ip  = (p  ) - (ip  ) as  defined  previously,  we  may  express  the 


n H 


n'h 


integral  equation  (6)  in  the  form  of  equation  (69) , where  the  potential  here  is 
defined  as 


P T(x*)  = - ip  (xA)  + 4>,(xa)  + 


Gq^dxdy  + 


fa°-^)Gnda  + p24 

h ' c *— 

i 


yds 


Gq°da  - F2  *p  Gq^vyds  . 
h ' c 


That  the  above  expression  for  the  potential  ij;  is  actually  identical  to  expression  (69a) 
may  be  verified  by  using  equations  ( 5 Id ) , (7)  and  (36a). 
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Figure  2:  The  correction  constant  k(B,6)  as  a function  of  the 
beam/length  ratio  B and  for  several  values  of  the 
draft/length  ratio  6. 
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